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ABSTRACT 

We present the luminosity function and color-redshift relation of a magnitude-limited sample of 145 mostly red 
field E/SO galaxies at z <J 1 from the DEEP Groth Strip Survey (GSS). Using nearby galaxy images as a training 
set, we develop a quantitative method to classify E/SO galaxies based on smoothness, symmetry, and bulge-to-total 
light ratio. Using this method, we identify 145 E/SOs at 16.5 < / < 22 within the GSS, for which 44 spectroscopic 
redshifts izspec) are available. Most of the galaxies with spectroscopic redshifts (86%) form a red envelope in 
the redshift-color diagram, consistent with predictions of spectral synthesis models in which the dominant stellar 
population is formed at redshifts z ^ 1.5. We use the tight correlation between V-I and Zspec for this red subset to 
estimate redshifts of the remaining E/SOs to an accuracy of ^^10%, with the exception of a small number (16%) 
of blue interlopers at low redshift that are quantitatively classified as E/SOs but are not contained within the red 
envelope. Constructing a luminosity function of the full sample of 145 E/SOs, we find that there is about 1.1-1.9 
magnitude brightening in rest-frame B band luminosity back to z ~ 0.8 from z = 0, consistent with other studies. 
Together with the red colors, this brightening is consistent with models in which the bulk of stars in red field E/SOs 
formed before z/or L5 and have been evolving rather quiescently with few large starbursts since then. Evolution 
in the number density of field E/SO galaxies is harder to measure, and uncertainties in the raw counts and their 
ratio to local samples might amount to as much as a factor of two. Within that uncertainty, the number density 
of red E/SOs to z — 0.8 seems relatively static, being comparable to or perhaps moderately less than that of local 
E/SOs depending on the assumed cosmology. A doubling of E/SO number density since z = 1 can be ruled out with 
high confidence (97%) if f2„, = 1. Taken together, our results are consistent with the hypothesis that the majority 
of luminous field E/SOs were already in place by z ~ 1, that the bulk of their stars were already fairly old, and that 
their number density has not changed by large amounts since then. 

Subject headings: cosmology: observations - galaxies: evolution - galaxies: luminosity function 



1. EMTRODUCTION 

In the local universe, E and SO galaxies are overwhelmingly 
found to be red galaxies that have smooth, centrally concen- 
trated surface brightness profiles and a high degree of elliptical 
symmetry. Giant E/SOs have surface brightness profiles domi- 
nated by the r'/"* law, which also characterizes luminous bulges 
of spiral galaxies (see Mihalas & Binney 1981 and references 
therein). 

The formation of E/SO galaxies is of interest because it is a 
sensitive barometer for the timescale of the formation of galax- 
ies and structure in the universe. There are two extreme the- 
ories for the formation of E/SOs. One is that E/SOs formed 
very early via monolithic collapse of the protogalactic gas at 
high redshift (the so-called monolithic collapse model; Eggen, 
Lynden-Bell, & Sandage 1962; Larson 1975). The other is that 
E/SOs formed via a continuous process of merging, as occurs in 
hierarchical models of structure formation in which the matter 
density is dominated by cold dark matter (CDM; Blumenthal 
et al. 1984; Baron & White 1987). For example, in a flat uni- 
verse with 0,„ = 1, some hierarchical merger models predict 
that as many as 50-70% of present-day ellipticals assembled 
later than z = 1 (Kauffmann, White, & Guiderdoni 1993; Baugh, 



Cole, & Frenk 1996a). However, this formation epoch is sen- 
sitive to the assumed values of cosmological parameters, with 
the major merger epoch occurring at higher redshift in universes 
with lower matter density (Kauffmann & Chariot 1998). Other 
physical processes, such as more frequent merging at higher 
redshifts, could also influence the prediction (e.g., Somerville, 
Primack, & Faber 2000). 

A key observational test of whether E/SO galaxies have 
evolved strongly at recent epochs is to measure the number den- 
sity and colors of carefully selected E/SO galaxies out to z ^ 1 . 
An L* early-type galaxy {Mb ~ -19.4 + 51og(/z); Marzke et al. 
1998) would have an apparent magnitude of / ~ 20-21 at z — 1 
depending on the details of its luminosity evolution. Currently, 
the determination of redshifts of / ~ 22 early-type galaxies is 
feasible with ground-based telescopes of aperture ;> 4 m (Koo 
et al. 1996; Lifly et al. 1995; Cowie et al. 1996; Cohen et al. 
1996, 2000). Moreover, previous studies show that morpholog- 
ical classification of faint galaxies is possible down to / ~ 22.5 
using HST images with exposure times ;>1 hour (Griffiths et 
al. 1994a, 1995b; Driver, Windhorst, & Griffiths 1995; Glaze- 
brook etal. 1995;Imetal. 1996, 1995a; Abraham et al. 1996; 
Brinchmann et al. 1998; Schade et al. 1999; Im, Griffiths, & 
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Ratnatunga 2000; Roche et al. 1996, 1998). Thus, in principle, 
a combination of HST imaging and spectroscopy using large 
ground-based telescopes can provide the data to constrain the 
number density evolution of E/SOs at z < 1. 

Testing for evolution in the comoving number density of any 
galaxy population requires selecting exactly the same kinds of 
objects at higher redshifts that comprise the target sample at the 
current epoch. In this paper, we aim to select galaxies that meet 
quantitative morphological criteria for being E/SOs. Galaxy 
colors redden rapidly after a starburst, approaching within 0.2 
mag of their asymptotic B-V values in less than 2 Gyr (e.g., 
Worthey 1994). This is also roughly the time scale for mor- 
phological peculiarities induced by starbursts to smooth out 
and disappear (Mihos 1995). Essentially, our selection criteria 
choose just those galaxies at each redshift that are morphologi- 
cally the most symmetric and smooth. We expect that galaxies 
selected this way have had the least star formation over the last 
few Gyr well past their last merger, and, hence, should occupy 
the red envelope of galaxies at each redshift. If galaxies reach 
the red envelope and stay there, our method measures how their 
numbers accumulate over time. 

Previous studies of the evolution of field E/SO number den- 
sities utilize various data sets and reach contradictory conclu- 
sions. To study distant early-type galaxies that did not have 
HST images, Kauffmann, Chariot, & White (KCW; 1996) use 
the data from the Canada France Redshift Survey (CFRS; Lilly 
et al. 1995) and select early-type galaxies from the red en- 
velope based on evolving model color curves in z vs. V -I. 
Based on {V /Vmax) statistics (Schmidt 1968), they claim a sig- 
nificant deficit of red early-type galaxies at z ~ 1 , with the num- 
ber density there being roughly one-third of the value at z = 0. 
They interpret this strong evolution as evidence for recent as- 
sembly and star formation in E/SO galaxies. However, subse- 
quent works have shown that the CFRS sample is deficient in 
red galaxies beyond z = 0-8, probably due to incompleteness in 
the redshift measurements, and that cutting the sample lower 
than that redshift removes all evidence for evolution (Totani & 
Yoshii 1998; Im and Casertano 2000). These works also note 
that using model color curves to select galaxies from the red 
envelope alone is risky; models by different authors vary, and 
the region identified as the red envelope depends heavily on 
the history of star formation assumed in the models. Possible 
systematic errors inV-I can also affect sample selection sig- 
nificantly. 

Both qualitative and quantitative classification methods have 
been applied to deep HST images to isolate samples of mor- 
phologically selected early-type galaxies. Schade et al. (1999) 
study CFRS galaxies that also have HST images. They choose 
elUpticals based on a good fit to an law together with an 
asymmetry parameter; no color cut is used. They see no drop- 
off in numbers out to z = 1, in apparent contradiction to KCW, 
but many of their most distant E/SOs are blue and lie below the 
color cut adopted by KCW — it is not clear that these objects 
would be classified as E/SOs if seen locally. A similar con- 
clusion regarding the lack of large number-density evolution of 
E/SOs is also reached by Im et al. (1999), who combine spec- 
troscopic redshifts available in the Uterature with existing HST 
images. E/SOs in Im et al. (1999) are selected as galaxies that 
have a significant bulge (B/T J> 0.3) and appear morpholog- 
ically featureless when visually classified. They present red- 
shift distributions of various galaxy types and find that those 
of E/SOs is consistent with no number density evolution out to 
z 1 (for similar works, see Roche et al. 1998; Driver et al. 



1998; Brinchmann et al. 1998). However, the number of early- 
type galaxies in both of these studies is small (about 40), and 
thus it is difficult to draw firm conclusions. 

Results from the number counts of larger sets of morpho- 
logically selected early-type galaxies are also consistent with 
little number-density evolution out to z = 1 (Im, Griffiths, & 
Ratnatunga 2000; Driver et al. 1996, 1998; Menanteau et al. 
1999). For example, Menanteau et al. (1999) study galaxies in 
deep HST archive exposures, supplemented by ground-based 
infrared Hh-K' near-infrared imaging. They classify galaxies 
morphologically using both quantitative (A/C method; Abra- 
ham et al. 1996) and quahtative (visual) methods. Lacking red- 
shifts, they model galaxy counts versus color and conclude that 
the number density of distant spheroids has not changed sub- 
stantially since z ^ 1 . However, this conclusion applies only 
if blue spheroids are included; a major deficit of red spheroids 
v/ithV -{H+K') > 2.0 is seen. Since models predict that these 
are precisely the colors of red spheroids at z ^ 1, this could 
imply a strong decline in the number of galaxies in the red en- 
velope at z > 1, although this does not put strong constraints 
on evolution at z < 1 (for similar works but different views, see 
Barger et al. 1999; Bershady, Lowenthal & Koo 1998; Broad- 
hurst & Bouwens 2000). However, modeling counts is noto- 
riously sensitive to the assumed luminosity function and star 
formation history, and essentially similar count data can also 
be fitted to models with substantial number density evolution 
(Baugh et al. 1996b; Im et al. 2000). Redshift information 
remains essential to break the degeneracy in the model predic- 
tions. 

Our method in this paper uses the comoving luminosity func- 
tion to characterize the number density of distant E/SO galax- 
ies. This method yields both a characteristic magnitude M* as 
well as the local number density normalization 0* . Variations 
in M* versus redshift are an additional measure of galaxy evo- 
lution; if ellipticals are just forming at z ^ 1, we would expect 
their luminosities to be very bright there (and their colors to be 
very blue). Luminosity evolution can also be tracked in other 
ways, for example, by studying zeropoint offsets in the size- 
luminosity relation or the fundamental plane, and color offsets 
can be measured using residuals from the color-magnitude re- 
lation. The absolute stellar ages of E/SO stellar populations can 
also be measured by fitting to broadband colors. Unlike the 
luminosity function method, these techniques all yield useful 
information with just a few objects, but their conclusions may 
in consequence be less general. 

As a group, such studies are converging on the consensus 
that distant E/SO galaxies show both luminosity and color evo- 
lution, and that these effects are are consistent with models in 
which the bulk of stars in E/SO galaxies formed before, in some 
cases well before, z ~ 1. Many studies of both field and cluster 
E/SOs indicate that the rest-frame B-band surface brightness of 
these galaxies is brighter at higher redshifts, consistent with the 
expected brightening of passively evolving stellar populations 
formed well previously (Schade et al. 1997, 1999; van Dokkum 
etal. 1998; Bender et al. 1998; Kelson et al. 1997; J0rgensen 
et al. 1999; Treu et al. 1999; Gebhardt et al. 2000). Clus- 
ter colors are basically consistent with this picture (e.g., Stan- 
ford, Dickinson, & Eisenhardt 1998), but some surveys find that 
roughly 30-50% of morphologically normal field E/SO galaxies 
are quite blue due to recent star formation (Schade et al. 1999; 
Menanteau et al. 1999; also see Abraham et al. 1999). Other 
field studies do not find many blue, bright early-type galaxies 
out to z ~ 1 (Im et al. 1996; Im, Griffiths, & Ratnatunga 2000; 
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Kodama et al. 1999; Franceschini 1998). The discrepancies be- 
tween these studies appear to be caused by differences in mor- 
phological classification and sample- selection criteria. Quanti- 
tative classification of E/SOs, such as we apply here, can help 
resolve these by offering an objective way of selecting E/SOs. 

The previous study of the E/SO luminosity function that most 
closely resembles ours is by Im et al. (1996). These authors 
use 376 visually classified field E/SOs at / < 22 from the HST 
Medium Deep Survey and other HST fields to construct a lu- 
minosity function (LP) to z 1 .2. They find that L* field E/SO 
galaxies brighten by 1-1 .5 magnitudes back to z ^ 1 , consistent 
with passive evolution models. They also find no significant 
number density evolution from z ~ 0.2 to z ~ 1. The sample 
used by Im et al. (1996) is the largest to date with complete red- 
shifts and Hubble-type classifications from HST images. How- 
ever, the redshifts are based on V - / colors with only a Umited 
number (23) of spectroscopic calibrators. The morphological 
classifications are also mostly based on qualitative visual meth- 
ods. 

The present paper follows the basic plan of Im et al. (1996) 
but with several important additions. Like that paper, we start 
with a sample of galaxies with HST images (in the Groth strip) 
and with spectroscopic redshifts {Zspec, from Keck) that are used 
to calibrate our photometric redshifts based on V -/. However, 
our sample of spectroscopic z's is much larger (44) and extends 
more uniformly to z ~ 1 . Second, we invest considerable effort 
in developing morphological classification criteria that isolate 
red-envelope galaxies to high accuracy. We show that the se- 
lection procedure works using the Keck redshift sample, which 
allows us to plot whether candidate galaxies actually lie on the 
red envelope; we tune the procedure so that they do. A final 
check by color cut in V -7 vs. 7 succeeds in identifying and 
correcting for a small fraction of blue interlopers. For the re- 
sultant red-envelope sample, V -7 is an excellent photometric 
redshift indicator, again demonstrated by the test sample with 
Keck redshifts. For the blue interlopers, the photometric red- 
shifts will be underestimated, but we find that they do not bias 
our results (< 10%) since their number is small, and they af- 
fect mostly the very faint end of the LF. For that reason, we do 
not try to eliminate the blue interlopers from our sample by the 
color criteria. This is useful for keeping our selection criteria 
simple. The net result is a technique that can identify E/SO can- 
didates that predominantly occupy the red envelope to a given 
magnitude limit using only V and I HST images. This is used to 
expand the final sample, which numbers 145 GSS galaxies, by 
including galaxies that lack spectroscopic redshifts. With this 
significant sample size, our study places a significant constraint 
on the number density and luminosity evolution of luminous 
field E/SO galaxies to z 1 . 

This paper is one of three papers on the evolution of early- 
type galaxies based on DEEP data. Results on the study of 
luminous bulges will be presented in Koo et al. (2000), and the 
fundamental plane of field E/SOs out to z ~ 1 will be explored 
by Gebhardt et al. (2000). Section 2 presents the basic observa- 
tional data for the present study. Section 3 develops the method 
for selecting E/SOs using HST images and explains how it is 
applied. Section 4 shows our fits for the luminosity function, 
and Section 5 provides a short discussion of the results. Con- 
clusions are given in Section 6. 

The Hubble constant, 77o, is quoted as /z = 7/o/(100km sec"' 
Mpc~'), and we adopt the value h = 0.7 throughout this paper. 

2. OBSERVATIONAL DATA 



Our goal is to identify moderate-to-high redshift red E/SOs 
and measure their redshifts photometrically from HST images 
alone. We proceed in two steps. First, we develop and tune 
our selection method for finding E/SO galaxies by testing it on 
a subsample of GSS galaxies having spectroscopic redshifts. 
This test allows us to verify simultaneously that we are se- 
lecting red-envelope galaxies and that our photometric redshift 
scheme works for such galaxies. Note that the completeness 
of this spectroscopic training sample does not matter since the 
completeness of the final sample is limited only by magnitude. 

The Deep Extragalactic Evolutionary Probe (DEEP) is a 
study of faint field galaxies using spectra obtained with the 
Keck 10-m telescope and high-resolution images obtained by 
the Hubble Space Telescope (HST). As a part of DEEP, we have 
extensively studied the Groth Strip, which is comprised of 28 
contiguous HST WFPC2 fields covering roughly 134 arcmin^ 
on the sky (Groth et al. 1994; Koo et al. 1996). In this section, 
we summarize the basic data from which our E/SO galaxies are 
selected. 

We adopt the Vega magnitude system as in Holtzman et al. 
(1995) and in previous DEEP papers. The zero-point of the 
magnitude system depends on the detailed shape of the filter 
and CCD response function. The 7-band magnitude measured 
through the F814W filter is simply referred to as "7". Like- 
wise, for the V-band magnitude measured through the F606W 
filter we will use the symbol "V". Readers are alerted that this 
y-band magnitude from F606W can be up to 1.0 magnitude 
brighter than the conventional Johnson V magnitude, depend- 
ing on the color and redshift of the galaxy. The 7-band mag- 
nitude with the F814W filter is roughly equal to the Cousins 
7 magnitude, with Ipsuw - Icousm — 0.04-0. 1 (Fukugita, Shi- 
masaku & Ichikawa 1995). Vega-calibrated magnitudes can be 
converted to the AB -magnitude system using the following re- 
lation: 7AB-7=0.44andyAB-y = 0.11 (Simardetal. 2000). 

2.1. Photometric data 

Galaxies were imaged by WFPC2 in both the F814W filter 
and the F606W filter, with exposure times of 4400 s and 2800 
s respectively (one deep field has total exposure times of 25200 
s in F814W and 24400 s in F606W). The limiting magnitude 
for the detection of objects is about 7 = 24, and about 3000 
galaxies are detected to this magnitude using FOCAS (Groth 
et al. 2000; Vogt et al. 2000). For the 7 and V magnitudes, 
we use the model-fit total magnitudes throughout the paper as 
described below. To obtain structural parameters and morpho- 
logical parameters, the surface brightness of each object is fit 
by a 2-dimensional bulge-ndisk model profile using the G1M2D 
software package of Marleau & Simard (1998) and Simard et 
al. (1999, 2000). Here, the disks are assumed to have an expo- 
nential profile, and the bulge profiles are assumed to follow the 
r4 law (de Vaucouleurs 1948). This procedure is similar to that 
of Ratnatunga, Griffiths, & Ostrander (1999) and Schade et al. 
(1995), and returns model-fit structural parameters that include 
total magnitudes (I,oi, V,ot), bulge-to-total light fraction (B/T), 
disk scale length (ro), bulge effective radius (r^.), half-light ra- 
dius (rfii), and position angles and ellipticities for both the bulge 
and disk components. Sizes are defined along the major axis of 
the model fit. Instead of listing each magnitude as Iiot or Vtott 
we omit the subscript "tot". For a more complete description 
of the basic photometric catalog of GSS galaxies, see Vogt et 
al. (2000) and Koo et al. (1996). For further description of the 
model-fit procedure and a complete fist of derived structural pa- 
rameters, see Simard et al. (2000). 
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2.2. Spectroscopic data 

Spectroscopic data were taken during the 1995-1998 Keck 
observing runs using multislit masks on the Low Resolution 
Imaging Spectrograph (Oke et al. 1995). Each mask contains 
about 30^0 objects. The blue part of the spectrum is covered 
with a 900-Une/mm grating, and the red part of the spectrum is 
covered with a 600-line/mm grating, giving us spectral cover- 
age from roughly 5000 to 8000 . Spectral resolution is FWHM 
=~ 2.5 and ~ 4 for the blue and the red settings respectively, 
providing sufficient resolution to measure internal kinematics 
of faint galaxies at moderate to high redshift down to a velocity 
width of (7 ~ 70 km sec"' . The high resolution also resolves 
the [O 11] doublet emission at 3727 , which enables redshifts 
to be measured from this single line alone. It also improves 
removal of the night sky lines, thus increasing the likelihood 
of identifying lines in the red part of the spectrum. For spec- 
troscopic observations, objects were mostly selected based on 
mean (V^/i^ +I^/i^)/2 magnitude (i.e., roughly an R magnitude) 
through a l."5 aperture, and aperture color iV-I)ap through a 
1 "0 aperture. A higher priority was given to the selection of 
red objects with (V-I)ap > 1 -75, which roughly corresponds to 
a passively evolving stellar population beyond z = 0.7: there- 
fore, our sampling strategy provides more spectra of red E/SOs 
at z > 0.7 than expected with a uniform sampling strategy. The 
typical exposure time for each mask was 1 hour in each wave- 
length setting, and objects that did not yield redshifts on the first 
try were attempted on subsequent exposures. The total number 
of redshifts is about 590, and the sample is 90% complete down 
to / = 23. This limit is about 1 .5 magnitudes fainter than that of 
4-m class surveys (e.g., CFRS) and comparable to other Keck 
telescope faint-galaxy surveys (Cohen et al. 1996; Cowie et al. 
1996). 

For further description of the spectroscopic sample, see 
Phillips et al. (2000) and Koo et al. (1996). 

3. SELECTION OF E/SO GALAXIES 

3.1. Overview 

The ultimate goal is to identify red, moderate-to-high red- 
shift E/SOs and to measure their redshifts photometrically from 
HST images alone. This section tests the adopted selection pro- 
cedure using simulated images of local galaxies and shows that 
the number of selected galaxies is reasonably stable versus red- 
shift. 

In detail, there are two selection criteria that are somewhat 
in conflict. On the one hand, we seek to restrict the sample 
to objects whose redshifts can be estimated from V-I alone. 
This requires that the selected objects lie on the red envelope 
of galaxies in color versus redshift, and thus that they be mor- 
phologically quite "pure." On the other hand, we seek a clas- 
sification system that yields stable numbers of galaxies as the 
quaUty of data declines at high z- This requires locating the data 
cuts such that equal numbers of objects are scattered out of the 
sample by errors as are scattered in. However, the in-scattered 
objects are in general not red and will not yield good photomet- 
ric redshifts. Being blue, they appear to have low reshift, and 
they inflate the number of E/SOs counted at low z. 

The first step in developing the selection procedure uses the 
spectroscopic redshift subsample to verify that the selection 
procedure finds mainly red-envelope galaxies, by plotting can- 
didate E/SOs with redshifts in V - / vs. z- Having finalized 
the selection procedure, we again use the spectroscopic red- 



shift subsample to calibrate the relation between spectroscopic 
redshifts and photometric redshifts determined from V-I color 
(for the red subset of galaxies). A final photometric sample is 
selected from the full database (lacking redshifts) and numbers 
and densities are calculated in Section 4. 

Traditionally, local E/SOs have been visually classified. 
However, visual classification is subjective, and boundaries are 
vague between morphological types (e.g., E/SOs vs. SO/a's). 
Naim et al. (1995) compared independent visual classifications 
of local galaxies by various experts and found that the average 
uncertainty over all types is 5T ^ 2, where T is the morpho- 
logical type defined by the RC3. The same uncertainty appUed 
to galaxies with T <Q (E/SOs). There are other difficulties as 
well; morphological features can be lost when spiral galaxies 
are viewed edge on, or when the number of resolution elements 
declines at high redshift and is insufficient to resolve spiral arms 
or other morphological details (so-called pixel smoothing). 

To lessen these difficulties, we establish a quantitative 
scheme that promotes a more objective and reproducible mor- 
phological classification. Such a scheme also enables us to use 
simulations to test efficiently the effects of pixel smoothing and 
low signal-to-noise (5/AO that afflict high-redshift data. 

The properties of E/SOs suggest that two parameters should 
suffice for the purpose of quantitative classification. One pa- 
rameter should describe how featureless and symmetric the 
appearance of the galaxy is, since E/SOs are characterized as 
smooth, featureless, symmetric objects. Another parameter 
should describe how the light is distributed within the galaxy, 
since E/SOs have a centrally concentrated surface brightness 
profile well fit by the r* law. 

As the two morphological selection parameters, we have 
therefore chosen to use bulge-to-total light ratio, B/T, to de- 
scribe the concentration of light, and the residual parameter, R, 
to describe the smoothness and symmetry of the galaxy mor- 
phology. Both parameters are byproducts of the 2-dimensional 
surface brightness profile fitting technique described in Section 
2.1. 

3.2. Bulge-to-total light ratio, B/T 

The quantity fi/r is the fraction of the Ught contained in the 
bulge component. Suppose that Lb and Lq are the total light in 
the bulge and disk respectively. Then, B/T is given as 



The quantity B/T is well correlated with the concentration pa- 
rameter (C), an alternative indicator of the concentration of light 
that has been used in previous studies of galaxy classification 
(e.g., Abraham et al. 1996). 

Kent (1985) shows that almost all SO galaxies in his 
magnitude- limited sample of nearby galaxies have B/T > 0.3. 
For E galaxies, Scorza et al. (1998) find that some of the 
28 E galaxies in their sample have B/T as small as ^ 0.5. 
Thus, adopting a B/T cut at around 0.4 is a reasonable choice 
for selecting E/SOs. Also, B/T is useful for excluding dwarf 
spheroidal or other late-type galaxies that may have a smooth 
appearance but have SB profiles closer to an exponential law 
(e.g., Im et al. 1995b and references therein). 

However, previous works have found that some late-type 
galaxies can have B/T > 0.4 and/or high C values (e.g., Kent 
1985). For example, about 20-30% of galaxies in the Kent 
(1985) sample with B/T > 0.4 are later than E/SOs. Therefore, 
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B/T (or C) alone is not sufficient to select E/SO galaxies, since 
early-type spirals and late-type galaxies with B/T > 0.4 would 
significantly contaminate the sample. For this reason we add a 
second parameter based on galaxy symmetry and smoothness. 

3.3. Residual parameter, R 

Our 2-dimensional surface-brightness fits provide a smooth 
and symmetric best-fit model for each galaxy. Subtracting the 
model image from the actual image yields a residual image. 
These residual images contain valuable information on galaxy 
morphology (Naim, Ratnatunga, & Griffiths 1997; Schade et al. 
1995). Figure 11] presents representative images of E, Sbc, and 
Im galaxies taken from the catalog of nearby galaxies of Frei 
et al. (1996). Corresponding residual images and 1-D surface 
brightness profiles along the major axis are also shown. The 
residual image of the elliptical galaxy (top) is smooth in com- 
parison to that of the spiral galaxy (middle), which shows spiral 
structure. The model fit to the Im galaxy (bottom) is poor, and 
many features can be seen in its residual image. 

To quantify the residual images, we use the residual param- 
eter, R, as defined by Schade et al. (1995, 1999). The quantity 
R is called the "asymmetry parameter" in these papers, but we 
adopt a different name, "residual parameter", since R measures 
how irregular the galaxy is, how prominent the spiral arms are, 
and how discrepant the galaxy surface brightness distribution is 
from the simple bulge-ndisk model, in addition to mere asym- 
metry. 

R is defined as 



R = Rj + R^^ 



with 



Ri 



Ra = 
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where Rij is the flux at pixel position (/, j) in the residual im- 
age, R]J^^ is the flux at (ij) in the residual image rotated by 
180 degrees, Bij and B]J^'' are similar quantities measured for 
background noise, and /,j is the flux at (ij) in the object im- 
age. The sum is over all pixels, and thus each pixel pair appears 
twice. We call Ra the "asymmetric residual parameter" since 
it is closely related to conventional asymmetry parameters, and 
we call Rt the "total residual parameter" since it measures the 
absolute strength of the residuals. The second terms in equa- 
tion ^ and (^ are approximate corrections for the contribution 
from random background noise (see Appendix A). When noise 
dominates the image, the first terms in Rj and Ra may become 
quite large since we may be adding up the absolute values of 
noise, and R may be overestimated. The noise is not negligible 
even when the S /N of the image is high (see also Conselice et 
al. 2000). Therefore, it is essential to include a background 
noise correction in equations (^ and (Q). 

An approximate expression for the error in R coming from 
the background noise is derived in Appendix A. The result is 



SR :^ 1.1(8/ N)- 



(5) 



This error does not contain contributions due to centroiding er- 
rors, which are present even at very high S/N. Conselice et 



al. (1999) have shown that the rms error in their asymmetry 
parameter. A, has a minimum value of 6A ^ 0.02 even with per- 
fect S/N. This directly translates to SR ~ 0.02. Therefore, for 
the final error in R we adopt the following relation; 



SR= J(0.02)^ + (l.7(S/N)-^)^. 



(6) 



In Schade et al. (1995, 1999), Rt and Ra are calculated 
within a 5 kpc radius from the center of each galaxy image. 
Calculating R within a fixed physical radius may not always be 
preferred since this would sample different parts of galaxies de- 
pending on how extended in physical size they are. To lessen 
such a problem, we calculate R within 2 r/,/. This has the ad- 
ditional advantage that we do not need to recalculate R when 
different cosmological parameters are assumed. 

As an illustration of the advantage of R over more conven- 
tional asymmetry parameters for selecting E/SO galaxies, con- 
sider a spiral galaxy with perfectly symmetrical spiral arms. 
Conventional asymmetry parameters, computed by subtracting 
the 180 degree-rotated image from the original image (Abra- 
ham et al. 1996; ConseUce et al. 2000; Wu 2000), would mea- 
sure the galaxy to be "symmetric" with no information what- 
soever about spiral arms. With R, one is able to quantify the 
prominent features caused by spiral arms in the residual image. 

3.4. Test on a local galaxy sample 

This section demonstrates that R and B/T can be used to se- 
lect E/SOs without substantially contaminating the sample with 
other galaxy types. Figure [l| shows B/T and R values below 
the residual image of each galaxy. The general trend is such 
that early-type galaxies have small R and large B/T, while late- 
type galaxies have large R and small B/T. 

Figure ^ plots R vs. B/T for the subset of 80 galaxies from 
Frei et al. (1996) that have enough background area necessary 
to determine a proper sky background subtraction. The Frei et 
al. sample contains a wide variety of Hubble types and is thus 
well suited for studying the relation between our quantitative 
morphological selection criteria and visual classifications. 

The red squares in Figure ^ represent galaxies with RC3 type 
less than or equal to -3 (E or E/SO), green triangles represent 
galaxies with T=-2 (SO), stars are for objects with -2 < T < 
(SO or SO/SOa), and the black crosses represent other types 
(T > 0). The box drawn in Figure || corresponds to the bor- 
der defined hy R < 0.08 and B/T > 0.4. Objects inside or on 
the border of the box are classified as QS-E/SOs ("QS" denotes 
"quantitatively selected"), and we find that nearly all of the se- 
lected objects are E/SOs with T <-0 (16 out of 17). Impor- 
tantly, there are almost no spiral or peculiar galaxies (T > 0) in 
the box (only one object). In fact, the residual parameter (R) 
cut by itself can provide a sample dominated by E/SOs. Thus, R 
alone can be used to define local E/SOs, but B/T will become 
important at high redshift when pixel smoothing degrades im- 
age quality (see Section 3.5). 

Note that roughly one-third of the Frei et al. galaxies with 
T <0 are missed with these selection criteria. The objects omit- 
ted are mostly of borderline type, with -2 < T < (SO, SO/SOa). 
Some of these could be included by loosening the selection cri- 
teria, at the expense of contaminating QS-E/SOs by non-E/SOs. 
We have examined images of the missed objects with T <Q and 
found that several have a non-smooth or irregular appearance, 
whereas the QS-E/SOs selected using the above criteria appear 
to be truly regular systems. Adoption of these stricter crite- 
ria is in keeping with our wish to minimize contamination by 
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blue interlopers. On balance, our scheme of selecting E/SOs is 
probably somewhat more conservative than the morphological 
typing of the RC3, and we seem likely to miss about 20-30% of 
E/SOs with RC type T < 0. We make use of this fact below to 
correct the counted numbers of E/SOs in local surveys. Effec- 
tively, our morphological cut roughly corresponds to T ~ -2, 
but we expect that more -2 <^T <^Q objects will be chosen if 
we apply our method on noisier images. 

3.5. Tests on "shrunken" local galaxy images 

We have shown in the previous section that our quantita- 
tive scheme is effective at selecting local E/SOs. However, one 
cannot blindly apply the above criteria to HST images of faint 
galaxies because of pixel smoothing. As we look at more dis- 
tant objects, each pixel of a given image samples a larger phys- 
ical area. Because of this effective smoothing, our/? (as well as 
other asymmetry parameters in the literature) is underestimated 
for galaxies with small apparent size. Note that B/T is much 
less susceptible to pixel smoothing since the bulge-ndisk fit pro- 
cedure incorporates the effects of pixel binning. However, it is 
not entirely free from systematic errors arising from either pixel 
smoothing or reduced S/N (see Simard et al. 2000). 

In order to see how important pixel smoothing and S/N are, 
we block-average images of the same 80 Frei et al. galaxies to 
simulate galaxy images with different half-light radii of roughly 
5, 3, 2, 1.5, and 1 pixel. Note that the apparent sizes of dis- 
tant galaxies in the GSS sample average rki,„,ed ^ 5 pixels at the 
magnitude limit of / ~ 22, and nearly all are larger than 3 pixels, 
so these tests are conservative (see Figure 6; also Simard et al. 
2000). Background noise is added to the resultant image so that 
the S/N of each simulated image is S/N ~ 40-80, comparable 
to GSS galaxies with / = 21-22 (see Figure Sample postage 
stamp images of simulated galaxies are available in Appendix 
B, along with their morphological parameters. In Figure]^, we 
show the input values of R (as derived from the simulated Frei 
et al. images with r/,/ = 5 pix) vs. output R for the simulated 
galaxies. Note that output r/,/ values from GIM2D are close to 
input values with a random error of order of Srui /rui ~ 0.08, and 
a slight systematic bias (~ 10% or more) for underestimating 
sizes when input sizes are <^ 2 pixels. What we actually have 
for GSS galaxies are "output" values, and we find that simu- 
lated Frei galaxies with input r/,/ ~ 1 .5 pixels have output r/,/ of 
rhLputput ^1.12 pixels. Similarly, for VhUnpu, ~ 1 pixels, we find 
rhi.putput ^ 0.72 pixels, and for r/,/,/„^„, ~ 2 pixels, rhi.outpui ^ 1 ■^ 
pixels. 

Importantly, output R does not change significantly from the 
input value (20% or less) when sizes of galaxies are sufficiently 
large (r/,/ > 2 pixels). When galaxy sizes become smaller than 
1.5-2 pixels, the global shift of R is significant (30%). 

Figure ^ is a similar comparison of output B/T vs. in- 
put B/T. When apparent sizes are very small (r/,/ ^1.5 pix- 
els), B/T values are again poorly determined. Simard et 
al. (2000) perform more extensive tests using artificially con- 
structed galaxies with various B/T values and find a similar 
result. However, only a small fraction of GSS galaxies in our 
sample have r/,; < 2 pixels; hence the effect of pixel smoothing 
on both the R and B/T cuts should be small. This is shown 
explicitly below. 

Figure ^ shows the R-B/T diagram for Frei et al. simulated 
galaxies with output r/,/ ~ 5, 3, 1.7, and 1.0 pixels. The selec- 
tion criteria shown on Fig. 5 (R < 0.08 and B/T > 0.4) select 
RC3 type E/SOs fairly well, although some E/SOs, especially 
with T > -2, are missed. However, the number of QS-E/SOs 



increases from 15 for r/,; = 5 pixels to 19 for r/,/ = 1 pixel, in- 
dicating that contamination from spiral galaxies becomes more 
important as galaxy size decreases. In order to compensate for 
this, we can try lowering the upper limit of the R cut for smaller 
galaxies. Adopting R = 0.07 for 2 < r/,/ < 3 pixels, R = 0.06 
for 1 < r/,/ < 2 pixels, and R = 0.05 for r/,; < 1 pixel is found 
to yield roughly stable numbers of E/SOs over the whole range 
of galaxy sizes in the simulated Frei et al. sample. However, 
we stress that the great majority (> 93 %) of GSS galaxies with 
/ < 22 have r/,/ > 2.5 pixels, as shown in Figure^, next section; 
thus any such reduction in R cut for small galaxies would not 
come into play for many objects. 

3.6. Final selection ofGroth Strip QS-E/SOs 

Using local galaxy images (Figure |l]), we have found that a 
constant boundary of R < 0.08 and B/T > 0.4 selects E and 
SO galaxies quite well without contaminating the sample with 
later galaxy types. However, when the object sizes are small 
(fhi < 2 pixels), pixel smoothing starts to wash away detailed 
morphological features, causing an underestimate of R values 
and a consequent overinclusion of galaxies. Also, as implied 
in equation (p|), errors in R increase cis S/N becomes smaller 
(S/N < 50). These effects are dealt with in this section. 

Figures ^ and ^ show r/,/ vs. / and S/N vs. / for galaxies 
in the GSS. Here, S/N is defined as the S/N within one r/,/ ra- 
dius of the object. The number of galaxies with r;,; < 3 pixels is 
not large but is not completely negligible; therefore we need to 
take into account the effect of pixel smoothing on R. Likewise, 
the median S/N approaches ~ 30 at our lower magnitude 
limit of / = 22. Using equation (5), we get an rms uncertainty 
SR ~ 0.06 at this brightness level, which is significantly larger 
than our rough estimate of the minimum scatter due to centroid- 
ing errors (6R ~ 0.02). Thus the effect of S/N decrease on R 
needs to be considered as well. After experimentation, we have 
adopted the final R cuts shown in Table 1 , which are a function 
of both r/,/ and magnitude. Tests below suggest that these R 
cuts, which are rather stringent, may be dropping some E/SOs 
at the faintest and smallest levels. However, we have retained 
them because they efficiently reduce the number of spurious 
blue interlopers while keeping the number of red E/SOs fairly 
close to intact. The efficacy of the adopted R cuts is examined 
next. 

E/SOs are the reddest galaxies in the local universe. If our se- 
lection method is good at identifying E/SOs and if these galaxies 
remain red at recent epochs (as will happen, for example, if the 
bulk of star formation occurs at z > 1), we expect our selected 
QS-E/SOs to populate a tight red envelope in the redshift-color 
diagram. Reassuringly, nearly all the objects selected this way 
are indeed the reddest galaxies at each redshift, as shown in 
Figure ||. This figure shows all 262 GSS galaxies with Zspec at 
16.5 < I < 22. Circle size is proportional to brightness, with 
the largest circles representing galaxies with / < 20, mid-sized 
circles objects with 20 < / < 21, and smallest circles objects 
with 21 < / < 22. Dashed lines indicate plausible ranges for the 
color of a passively evolving stellar population formed at very 
early times. The upper dashed line represents a 0.1 -Gyr burst 
model with 2.5 times solar metallicity, a Salpeter IMF (0.1 to 
135 Mq), and z/or =11. The lower dashed line represents the 
same 0. 1-Gyr model but with 0.4 times solar metallicity. To al- 
low for color errors, we have added or subtracted 0.15 mag in 
V-I to and from the upper and lower dashed lines. 

Panel b) of Figure pTshows the 84 galaxies in the spectro- 
scopic training sample that satisfy the R cut. The great majority 
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fall within the plausible color range of the passively evolving 
stellar population (51 out of 84 galaxies). By the same token, 
only 33 out of 186 galaxies outside the red color boundaries 
in panel a) are selected as low-/? galaxies. Figure ^ likewise 
shows the 77 galaxies that satisfy B/T > 0.4. Most of these 
again turn out to lie within the red color boundaries (50 out of 
77 objects), and only a small fraction of blue objects below the 
red boundary have high B/T (27 out of 186 galaxies). 

Finally, Figure ^ presents the 44 galaxies that satisfy both 
the R and B/T cuts. In addition to the previous dashed lines, 
we also plot a solar-metallicity model with three different for- 
mation redshifts (11, 2, and 1.5). Now, only 6 of 44 selected 
objects (15%) lie outside the red color boundaries; these are the 
"blue interlopers" referred to previously (filled circles). The re- 
mainder of the sample is found to follow the expected redshift- 
{V -I) tracks of passively evolving stars. The nature of the blue 
interlopers is intriguing. Preliminary analysis of their spectra 
(Im et al., in preparation) shows that most have strong, nar- 
row emission lines, suggesting that they are low-mass starbursts 
rather than massive star-forming E/SOs (Im et al. 2000, in 
preparation); they may be similar to the Compact Narrow Emis- 
sion Line Galaxies (CNELGs; Koo et al. 1994; Guzman et al. 
1997; Phillips et al. 1997). 

We now vary the selection rules to see how the final sam- 
ple depends on the precise criteria used. Figure |^ shows V — I 
vs. redshift for samples selected based on R and B/T cuts that 
are slightly different from those adopted in Figure ^. The two 
figures on the left (panels a and c) show the results of loosen- 
ing the cuts. The number of blue interlopers significantly in- 
creases, from 6 in Figure ^ to 11-12 here, while the number 
of selected red QS-E/SOs increases by only 1^. The two pan- 
els on the right (b and d) show the effect of stricter cuts. The 
number of blue interlopers is not significantly reduced (4 here 
vs. 6 formerly), while the number of desirable red QS-E/SOs is 
decreased significantly, by 7-10 objects. Thus, the original cuts 
seem about optimal. 

Interestingly, there may be a bimodality in the colors of GSS 
galaxies such that the color distribution at a given redshift is 
double-humped. A hint of this is seen in Figure || and has been 
remarked on previously (Koo et al. 1996) (but note that no such 
feature is seen in CFRS data; Lilly et al. 1995). Such a hump 
would clearly help in the selection of red E/SOs; setting crite- 
ria to cut in the valley would mean that object selection would 
be less sensitive to slight changes in the selection criteria. This 
approach can be tried in future if color-redshift surveys confirm 
the presence of the double humps. 

We next discuss likely systematic errors in the counted num- 
bers of E/SO galaxies at high redshift. Figure 1^ shows images 
of all selected GSS QS-E/SOs with Zspec^ ordered by redshift. 
Blue interlopers (objects lying outside the red bands of Figure 
H) are separately presented at the end of the sequence. The 
pixel values of the galaxy images are roughly square-rooted 
(more exactly, rescaled by the jV^h power of their values, as 
used in Frei et al. (1996)). We find this scaling to be effec- 
tive in bringing up faint details at low surface brightness, while 
making the bulge component look reasonably distinct when 
B/T > 0.4. However, for eyes accustomed to looking at the 
linearly scaled images of most astronomical atlases, the scaling 
used here might make B/T ^ 0.5 objects appear rather disk- 
dominated. To avoid this confusion, we add similarly scaled 
images of local E/SO galaxies from Frei et al. (1996) in the two 
bottom rows of Figure |l^. The visual comparison between the 
local E/SOs and the Groth Strip QS-E/SOs confirms that the lat- 



ter truly resemble the appearance of local E/SOs. Thus, aside 
from the blue interlopers (15%), which will all appear at low 
zs from their photometric redshifts, the present classification 
scheme admits at most few additional spurious spirals and pe- 
culiar galaxies and is thus not likely to overestimate the number 
of distant E/SOs by even a small percentage. 

For the reverse comparison. Figure 111] shows images of 64 
galaxies with Zspec that lie within or close to the red color 
boundaries of Figure ^but that do not meet the R or B/T cuts. 
Such objects could be real E/SOs that are improperly being lost. 
Comparison of this figure with Figure |l^ shows that red non- 
selected E/SOs actually have a much higher frequency of non- 
smooth morphological features (e.g., spiral arms, asymmetric 
nuclei), which are not apparent in selected QS-E/SOs. Many 
also turn out to be edge-on galaxies with low R but also with 
low B/T. However, a significant number of the non-selected 
galaxies are indistinguishable visually from the QS-E/SOs of 
Figure |lO[ We find about 10 such objects in Figure [iT], of which 
9 lie beyond z = 0.6. If these are truly E/SOs, they should be 
added to the 24 QS-E/SOs in that redshift range from Figure |l^ 
which would mean that our numbers of high-redshift E/SOs are 
^ 35% too low. These objects might overlap at least in part 
with the borderline SO-SO/a's missed in the test of local objects 
using the Frei et al. catalog in Figure ^. 

Clearly, resolution and S/N effects can work both for or 
against selecting E/SOs, but Figure [ll] suggests that, in our 
method, they seem to work against picking E/SO galaxies at 
z > 0.6 but do not seem to affect the low-redshift E/SO selec- 
tion very much. Only a few of z < 0.6 red non-E/SO galax- 
ies in Figure 111] would resemble QS-E/SOs at z > 0.6 if their 
V-band images were reprocessed to appear like z > 0.6 galax- 
ies. The great majority 50 %) of red non-E/SO galaxies at 
z < 0.6 are edge-on galaxies with negligible bulge component 
(B/T <, 0.2), while only - 25% of red, non-E/SO galaxies at 
z > 0.6 are in such category. This supports the above idea, and 
also suggests that red galaxies at lower redshifts can be more 
easily contaminated by dust-extinguished edge-on disks than 
red galaxies at higher redshifts. Again, the point is to establish 
that our counts are not likely to M«t/erestimate the number of 
distant E/SO galaxies by more than ^ 35%. 

Finally, we note that our measured values of R and B/T are 
derived from the observed /-band image, which, for our sample 
at z 0.8, corresponds to a rest-frame B band. Since our local 
galaxy comparison sample is observed in the B band, the mor- 
phological K-correction should be minimal when 0.6 < z < 1 .2. 
For galaxies at z < 0.6, however, this could bias object selec- 
tion because R and B/T will be estimated at redder rest-frame 
wavelengths. Since bulges are redder than disks and local- 
ized star formation is less prominent at redder rest-frame wave- 
lengths, the expectation is that B/T would be overestimated 
and R would be underestimated when measured in rest-frame 
V rather than B, and that consequently more objects would be 
selected as E/SOs at z < 0.5. To test this, we have compared 
B/T and R values measured in V vs. / and find there is no 
strong difference as long as both V and / sample light above 
rest-frame 4000 . Re-selection of the sample at z < 0.5 using 
V-band images rather than /-band shows further that a V-band 
selected sample would be almost identical to the /-band E/SO 
sample. We have attempted to estimate the rest-frame B-band 
B/T hy applying K-corrections from Gronwall & Koo (1996) 
and Gebhardt et al. (2000) separately to bulge and disk com- 
ponents, and confirm the above claim that rest-frame B-band 
B/T is nearly identical whether estimated from observed V or 
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/. Therefore, we believe that morphological /T-correction is not 
an important issue here. 

3.7. Selection ofGroth Strip QS-E/SOs without spectroscopic 

redshifts 

Our spectroscopic observations do not cover the entire Groth 
Strip, so we can more than triple the sample size by estimating 
redshifts photometrically for galaxies in regions where spectro- 
scopic data are not available. For the GSS as a whole, the only 
photometric information we can use are /, V, and (V-/). Due 
to the wide range of color space spanned by various types of 
galaxies, it is not feasible to estimate redshifts for all types of 
galaxies in the GSS using this limited photometric information. 
Nevertheless, it is possible to get reliable redshifts with {V -I) 
only if we focus on a sample of E/SOs preselected by morphol- 
ogy, since the previous analysis has shown that, for them, iV-I) 
color and redshift are very well correlated (cf. Figure The 
correlation is virtually perfect at / < 21, where the R and B/T 
cuts select E/SOs with a tight Zspec-{y -I) relation. However, at 
/ > 21, blue interlopers make photometric estimate of redshifts 
more challenging. In order to exploit the tight Zspec-iY -I) rela- 
tion for accurate photometric redshifts, we exclude blue inter- 
lopers from consideration and use the remaining QS-E/SOs to 
fit Zspec with V -I and / polynomials. We obtain the following 
relation: 

Zphot=Ax+A2l+A^f+A4f+AsIy(.{V-r) + 

+A(,{y-r)+AAV-if+A^(v-i)\ (7) 

where the coefficients are Ai = -2.7872 x 10"^ A2 = -1 .7700 x 
10-1, A3 = 1.6784 X 10-2, A4 = 4.3770 x A5 = 2.7295 x 
10-2, Afi = 2.0530 X 10-^ A7 = 2.1158 x 10-^, and Ag = 1 .6606 x 
IQ-^. The quantity Zphot estimated with equation (Q) appears 
to underestimate systematically the true redshift by a small 
amount at z ^ 0.8. For that reason we make the following small 
correction when Zpiwt > 0.8: 

Zcorrected ~ 

1.333(z-0.8) + 0.8. (8) 

Including terms in / as well as V- / in equation (|^) reduces the 
residuals by about 10%. 

Figure [l^ compares Zspec vs. Zphot for the GSS QS-E/SOs 
having spectroscopic redshifts (blue interlopers excluded). The 
RMS of Zphot vs. Zspec is about 10% but increases at the high- 
est redshifts, as in Figure |l3|. The lines there indicate the rms 
of Az = Zspec -Zphot VS. Zspec, and wc adopt this as the error 
of Zphot- This error envelope will be used later in the estimate 
of luminosity function parameters and in tests with Monte- 
Carlo simulations for checking Malmquist-like bias. The er- 
ror increases beyond z > 0.8 due to the fact that the main Zphot 
indicator — the 4000 break in the continuum of the spectral en- 
ergy distribution — passes through the F814W passband. How- 
ever, the combination of the following two facts makes V -I 
and / together useful for estimating redshifts to reasonable ac- 
curacy (< 15%) even at z > 0.8. First, the observed z vs. V -I 
relation is not completely flat beyond z > 0.8, contrary to the 
predictions of the passive evolution models plotted in Figure ^. 
The color-magnitude relation is at least partly responsible for 
this — the magnitude and redshift limits we adopt make only the 
intrinsically brightest, and thus the reddest, objects detectable. 
This effect acts to increase the average color vs. redshift, even 
when the color of any given galaxy would remain flat. Second is 



the familiar fact that, at fixed intrinsic V — I, dimmer-appearing 
galaxies are farther away. Thus apparent magnitude is by itself 
an indicator of redshift, independent of color The fit at z > 0.8 
is based on more than 15 E/SOs with known Zspec in this range; 
therefore, our Zphot can be considered reliable within the esti- 
mated errors even at 0.8 <z< 1.2. 

At low redshift (z < 0. 1), there is a second concern that small 
errors in photometric redshift (e.g., 6z = 0.05) lead to large er- 
rors in absolute magnitude. The errors shown in Figure |l3| im- 
ply that redshift errors remain fractionally small (10%) even at 
very low redshift. However, in practice the errors are poorly 
known below z <^ 0.3 since there are only two galaxies in this 
redshift range. To check for a potential bias due to the effect of 
low-redshift errors, we repeat the LF analysis below, increasing 
the lower redshift cut to z = 0.2, and show that this has little 
effect. 

A cautionary remark must be made regarding the photomet- 
ric redshifts of blue interlopers. Since the number of blue in- 
terlopers is small (<15% of QS-E/SOs), we do not try to ex- 
clude them from the sample using additional color cuts. How- 
ever, redshifts for the blue interlopers are underestimated using 
equation (^. Fortunately, with these redshifts, blue interlopers 
tend to be the faintest QS-E/SOs at a given V - / color, and thus 
they influence only the faintest part of the luminosity function 
at low redshift. Figure |l^ which plots / vs. V -I for various 
samples, sheds further insight into the number of blue interlop- 
ers. The squares refer to the Zspec sample. Thick squares show 
the red QS-E/SOs, while thin squares indicate blue interlop- 
ers as defined previously in Figure The lines represent the 
color-magnitude relation for a passively evolving elliptical with 
Mb = -18.3 and Mb =-20. 14 (L*), assuming ^0 = 0.1, /; = 0.7, 
solar metallicity, S alpeter IMF, and z for = 5 . The maj ority of red 
QS-E/SOs lie above this line, while all but one blue interloper 
in the Zspec sample lie below the line. The circled points denote 
the additional galaxies in the Zphot sample. As no independent 
redshifts are available for them, we do not have firm knowledge 
of which ones are blue interlopers. However, the circles lying 
below the line are candidate blue interlopers according to Zphot', 
there are 14 of these, among 101 objects, similar to the 6 inter- 
lopers out of 44 objects in the Zspec sample. Redshifts of both 
kinds of interlopers are likely to be severely underestimated, 
and we find that all of them have Zphot < 0.42 and Mb > -18.2. 
Thus, the blue interlopers probably overestimate the faint end 
of the LF at low redshift. 

Figure [l^ shows images of 98 out of these 101 QS-E/SOs 
in the range 16.5 < / < 22 in the Zphot sample, including blue 
interloper candidates. All red QS-E/SOs are presented, and 1 1 
out of the 14 candidate blue interlopers are shown at the end 
of the figure. We can again inspect the images of these objects 
as a sanity check for spurious late-type galaxies and find that 
contamination by late-types and peculiars is very small; by eye, 
only 2 out of 98 galaxies look mis-selected. Thus, aside from 
blue interlopers, the likely overestimate of distant E/SO galax- 
ies is again very small, even in the Zphot sample. (The opposite 
test of looking for objects missed among "red" galaxies, which 
we performed for the Zspec sample, is impossible here because it 
requires a spectroscopic redshift to define a "red" galaxy.) The 
ellipticity distribution of the GSS QS-E/SOs is presented in the 
next section, which further shows that they are similar to local 
E/SOs. With the addition of the Zzphot sample, we have a final 
sample of 145 QS-E/SOs at 16.5 < / < 22.0. Information on 
these QS-E/SOs is fisted in Table 2. 

Figures |l6| and |l7| show redshift vs. /,o, and redshift vs. Mg 
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(K-corrected only) for all GSS galaxies. QS-E/SOs are plotted 
as squares, and blue interlopers are plotted with triangles. Thick 
symbols represent the spectroscopic sample, and thin symbols 
represent the photometric redshift sample. Small crosses in 
Figure [l^ are the remaining galaxies (non-E/SOs) with spec- 
troscopic redshifts. Also plotted in Figure [l^ are lines of three 
different values of constant Mb, with (solid line) and without 
(dashed line) luminosity evolution. Note that Mb = -20. 14 rep- 
resents the L* of local E/SOs according to Marzke et al. (1998). 
The parameters for the open universe are adopted, and for lumi- 
nosity evolution we assume E{z) = 1 .7 x z,a s derived from our 
LF analysis for the open universe (see next section). In Figure 
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we plot only lines for Mb = -20. 14. 
A striking feature in Figures |l6| and |l7| is that there seem to 
be too many L> L* E/SOs beyond z > 0.6 if the no-evolution 
line is used as a reference. This overabundance is not observed 
when we count the number of L > L* E/SOs with respect to the 
evoliving-luminosity line, and this can be considered a qualita- 
tive indicator of luminosity evolution. In the analysis of the LF 
below, we will quantify the amount of luminosity and number 
density evolution in detail. A second important feature is the 
apparent lack of L > L* E/SOs at z < 0.2, but this can be at- 
tributed to the bright magnitude limit we adopted (/ = 16.5, see 
Figure |l6|). In the LF analysis below, we adopt a default red- 
shift range of 0.05 < z < 1.2, despite the fact that few L> L* 
galaxies are seen at z < 0.2. The techniques we use are adaptive 
enough to adjust for this, but just to check, we try increasing 
the lower redshift cutoff and confirm that there is little effect. A 
more serious deficiency of galaxies might also exist at z > 1.0, 
but one that cannot be explained simply by magnitude limits. 
We will again vary the upper redshift cutoff and find that our 
results are slightly more sensitive to this upper cut. 

3.8. Ellipticity distribution 

The ellipticity distribution of elliptical galaxies is known to 
be quite different from that of spiral and SO galaxies, and thus 
can be used as yet another independent check on the Hubble 
types of the selected sample. Most ellipticals look round or 
football-shaped with modest ellipticities, or equivalently, large 
axis ratios. On the other hand, the ellipticity distribution of SO 
galaxies is peaked at e ^ 0.5, while that of spirals is nearly flat at 
almost all ellipticities (Sandage, Freeman, & Stokes 1970; van 
den Bergh 1990; Fasano & Vio 1991; Franx, Illingworth, & de 
Zeeuw 1991; Lambas, Maddox, & Loveday 1992; J0rgensen & 
Franx 1994; Andreon et al. 1996; Dressier et al. 1997). 

Figure |l8], presents the ellipticity distribution of the 145 GSS 
QS-E/SOs with 16.5 < / < 22 as the thick histogram {Zphot sam- 
ple included). Ellipticities are measured at the 2 isophote; 
they typically increase slightly with radius, but the increase is 
not large beyond r ^2 r^. Also plotted are the ellipticity distri- 
butions of local Es (dashed line) and SOs (dotted line) in nearby 
clusters, taken from Dressier (1980). The thin line shows the 
combined nearby E and SO ellipticity distributions for a model 
with a relative SO fraction of 40%. As a reference, we also 
plot an ellipticity distribution of spirals with a dot-dashed line 
(Lambas et al. 1992). 

Since the ellipticity distributions of Es and SOs are distinc- 
tive, it is worthwhile trying to separate the two types of galaxies 
in the GSS sample. Since we do not have a reliable scheme to 
distinguish Es from SOs quantitatively, we instead divide the 
sample above and below Mb < -20. According to the LFs 
of local field E/SOs (Mai'zke et al. 1994) and cluster E/SOs 
(e.g.. Dressier et al. 1980), Es ai^e more abundant than SOs at 



Mb -20, while SOs are more abundant than E's at Mb <; -20. 

Figures 19 and ^ show the ellipticity distribution of Groth 
Strip QS-E/SOs divided this way. To estimate the absolute 
magnitude of each object, we assume an open universe with 
n„, = 0.2 and h = 0.7; passive luminosity evolution (as derived 
from the luminosity function in the next section) and the K- 
correction are both taken into account. The resultant elliptic- 
ity distribution of luminous Groth Strip QS-E/SOs is well fit- 
ted with a model distribution dominated by Es, amounting to 
^ 80% of the sample, while the ellipticity distribution of the 
faint sample resembles a model distribution dominated by SOs, 
amounting to ^ 70%. Thus, the combined ellipticity distribu- 
tion of local Es and SOs reproduces that of the Groth Strip QS- 
E/SOs fairly well, and neither the ellipticity distribution of local 
Es alone nor that of SOs alone is a good fit. Moreover, neither 
bright nor faint Groth Strip QS-E/SOs are consistent with the 
ellipticity distribution of local spirals; late-type galaxies do not 
therefore appear to contaminate our sample significantly. 

3.9. Comparison with previous studies 

Studies of early-type galaxies by Brinchmann et al. (1998) 
and Schade et al. (1999) have included objects from a part of 
the GSS. Using these overlaps, we compare our selection crite- 
ria with theirs. Brinchmann et al. (1998) use the AC-system, 
which originates from Abraham et al. (1996). The AC-system 
classifies galaxies based upon their asymmetry (A) and con- 
centration parameter (C). According to this scheme, early-type 
galaxies are selected largely based on the concentration param- 
eter, C, which correlates well with B/T. As expected from 
Figure 5, we find that the AC method tends to include some 
later-type galaxies with large C (or B/T). Specifically, of the 
14 AC-classified early-type galaxies in the Groth Strip, 9 are 
classified as QS-E/SOs by us, while 5 of them are classified as 
non-QS-E/SOs (for example, 074_2237 and 073_3539 in Figure 
11). In contrast, none of our Groth Strip QS-E/SOs are classified 
as non-E/SOs by Brinchmann et al. (1998). Thus, we conclude 
that our method is somewhat more conservative than the AC 
method in picking up only morphologically featureless E/SOs. 

Schade et al. (1999) compiled a list of elliptical galaxies 
at z < 1. Their criteria are that the galaxy should have a 
law-dominated profile and R < 0.1, slightly looser than the 
R <; 0.08 criterion adopted here. Eight ellipticals from Schade 
et al. (1999) are found in the our sample. Of these, 5 are clas- 
sified by us as QS-E/SOs, while one more is a blue interloper. 
This shows again that our object selection criteria are more con- 
servative than those of Schade et al. (1999), as expected from 
the tighter R cut. The Schade et al. sample also has a much 
larger scatter in {V -I) vs. z than ours, probably due in part 
to the looser R cut and possibly also to larger errors in their 
ground-based photometry. 

3.10. Selection errors and biases for the distant sample 

This section summarizes previous discussions of the selec- 
tion effects and adds some new tests to produce an overall esti- 
mate of count uncertainties. Here we consider internal errors in 
our own counts only — errors in matching to local E/SO counts 
are considered in Section 3. 

Here is a summary of the tests: 1) Varying R and B/T thresh- 
olds within plausible limits modulates the absolute number of 
counted galaxies by ±25%. 2) Varying the thresholds as a func- 
tion of galaxy size affects the counts of small galaxies. The fi- 
nally adopted, conservative thresholds lose about 20% of Frei et 
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al. simulated E/SOs at the smallest radii (<^2 px), which affects 
low-luminosity and distant galaxies the most. 3) Inspecting the 
Zspec sample visually for all conceivable interlopers reveals no 
new objects other than the known blue interlopers (6 out of 20 
objects below z < 0.6), and the Zphoi sample gives essentially 
identical results. Thus, the counts beyond z = 0.6 are unlikely 
to be biased too high by inclusion of late-type or peculiar ob- 
jects, while the nearby counts below z = 0.6 may be roughly 
30% too high due to the inclusion of blue interlopers in low- 
luminosity bins. 4) Inspecting the Zspec sample visually for all 
conceivable omitted E/SOs reveals only one or two new galaxies 
that might be added to the 14 existing galaxies at low redshift 
{z < 0.6), but turns up an additional 9 objects that might plau- 
sibly be added to the 24 existing objects at z > 0.6. Thus, the 
counts beyond z = 0.6 may be biased too low by about 30% due 
to omission of valid E/SO galaxies, while the nearby counts do 
not appear to be missing such candidates. 

The overall conclusion from these tests is that the counts are 
uncertain at the 30% level; they may be biased a little too high 
by this amount at faint absolute magnitudes for nearby redshifts 
(z < 0.6), and a little too low by a similar amount at all absolute 
magnitudes for distant redshifts (z > 0.6). 

4. VOLUME DENSITY AND LUMINOSITY EVOLUTION OF GROTH 
STRIP QS-E/SO'S 

In this section, we construct the luminosity function (LF) of 
distant Groth Strip QS-E/SOs and derive constraints on their 
number density and luminosity evolution. The entire sample 
including Zphot is used unless otherwise noted. 

4. 1 . Luminosity function 

To derive the evolution in luminosity and number density 
from the LF, we adopt two different approaches. One is to de- 
rive the LF parameters of the sample divided into two different 
redshift intervals (0.05 < z < 0.6, and 0.6 < z < 1.2); evolu- 
tion is then measured by comparing low- and high-redshift LFs 
to one another and to the local LF. A second approach intro- 
duces parameters for evolution and fits the (evolving) LF for 
the whole sample simultaneously. For the first method (Method 
1), we use the l/VmoA technique (Schmidt 1968; Felten 1976; 
Huchra & Sargent 1973; Lilly et al. 1995a) to estimate the LF 
in magnitude bins at low and high redshift, and then apply the 
method of Sandage, Tammann, and Yahil (STY, 1979) to esti- 
mate the LF parameters (with normalization provided by yet a 
third method). Method 1 is identical to the approach adopted 
by Im et al. (1996). For the second method (Method 2), we 
follow an approach similar to that of Lin et al. (1999), in which 
luminosity and number density evolution are each parameter- 
ized versus redshift, and these parameters are then solved for 
together with other LF parameters using the whole sample si- 
multaneously. 

4.2. Method 1 

4.2.1. 1 /V,nax method 

In the 1 /Vmax method, each galaxy in the sample is assigned a 
1 /Vnifl.i value, where Vmax is the maximum volume within which 
the galaxy would be observable under all relevant observational 
constraints including magnitude and redshift limits. The quan- 
tity V„,ax is calculated as 

/■minizi.Zml) 

ymax= / {dV/dz)dz, (9) 



where z\ and zi are the lower and upper limits of the redshift 
interval for which the LF is being calculated, nii and mi are 
the apparent magnitude limits of the survey, Zm\ and Zmi are 
redshifts where the galaxy would be located if it had apparent 
magnitude m\ and m2 respectively, and dV /dz is the comoving 
volume element per unit redshift interval. 

An absolute magnitude of each galaxy in the F814W pass- 
band {Mi) is calculated as 

Ml = m,-5\og,^,{ddz))-25-Ki{z), (10) 

where diiz) is the luminosity distance in Mpc and ;«/ is the 
apparent magnitude of the galaxy. For the /T-correction, we 
use the present-day model SED which was used in the color 
cut in Section 3.6 (i.e., the BC96 model with 0.1 Gyr burst, 
Zfor = 5, solar metallicity, and Salpeter IMF). This /T-correction 
is very similar to the A'-correction used for Es in Gronwall & 
Koo (1995), the standard set of /T-corrections used in previous 
DEEP publications. The difference between the two is roughly 
0.08 mag at z > 0.5 and ~ 0.3 mag at z ^ 0.2, with the adopted 
K-correction underestimating the luminosity in both cases. To 
obtain the rest-frame Mb magnitude, we add 2.17 mag to M/ 
since the (B-I) color of the model E/SO SED at z = is 2.17. 

Our sample of QS-E/SOs is magnitude-limited at 16.5 < / < 
22, as described in the previous section. Since the number of 
E/SOs with z <; 1 is small with a faint apparent magnitude limit 
of / = 22 (See Figure |l^, and since the accuracy of Zphot rapidly 
drops due toV-I color degeneracy beyond z^ 1 , we restrict the 
redshift interval to 0.05 < z < 1.2. With these selection criteria, 
the total number of QS-E/SOs is 145; the number with Zspec is 
44. When there is luminosity evolution, the real /Tz-correction 
should include the luminosity dimming term, E{z). This would 
change z,n2 and Zmi, thus affecting V„,ax as derived from equa- 
tion (ph. Intrinsically bright galaxy samples are nearly volume- 
Umited (i.e., miniziyZmi) = Zi and max(zuZmi) = Zi), except at 
very low redshift (z <^ 0.2). However, the volume at z < 0.2 is 
small compared to the remaining volume (e.g., 0.2 < z < 0.6), 
and the evolutionary correction itself is small at low redshift; 
thus, the evolutionary correction does not significantly affect 
the bright end of the LF (< 0.05 dex in density). As a result, 
1 /Vmax tends to become bigger with negative E{z). However, 
the level of change in 1 /Vmax is only ~0. 1 dex. 

Galaxies are then divided into different absolute magnitude 
bins, and the LF value for the jth bin is calculated as the sum 
of all 1 /Vmax values of galaxies belonging to that bin, i.e., 

0(M), fl'M=S (11) 

y/nax 

4.2.2. STY method 

To estimate the parameters of the LF, we use the STY 
method (Sandage, Tammann & Yahil 1979; Loveday et al. 
1992; Marzke et al. 1994, 1998; Efstathiou, Ellis, & Peter- 
son 1988; Willmer 1997), assuming that the LF is described by 
the Schechter form (Schechter 1976): 

(/)(M)(/M = 0.921 0*ji:"-''exp-^t/M (12) 

where ;c=10"-4('^*-'*^>. 

The Schechter function has three free parameters ((/>*, M*, 
and a); (j)* is for the density normalization, M* indicates the 
characteristic luminosity of the distribution where the number 
density of bright galaxies starts to fall off, and a is the slope of 
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the faint end of the luminosity function. For the LF of all nearby 
galaxies, these parameters are estimated to be = 0.01-0.02 
Mpc-^ M* = -19.0 to -19.7 + 51og(/!) B mag, and a = -0.8 to 
-1.4 (Marzke et al. 1998; Lin et al. 1996; Zucca et al. 1997) 

In the STY method, the luminosity function parameters are 
found by maximizing the probability of the observed data, and 
hence by maximizing the following likelihood function: 



InL = Sin Pi, 



(13) 



where is the normalized probability of finding galaxy / with 
absolute magnitude M, at redshift z in a magnitude-limited sur- 
vey. The normalized probability is given by 



Pi = mdi 



■miniMinaxMi) 



(j){M)dM. 



(14) 



max[M,„i„ ,M\ ] 



Here, M, is the absolute magnitude of the object, M\ and M2 
are the brightest and faintest absolute magnitude limits of the 
sample, and M,„ax and M„„„ are the maximum and minimum 
absolute magnitudes observable at redshift z given the apparent 
magnitude limits of the survey. In our analysis, we do not re- 
strict Ml and M2 so these quantities are irrelevant here. Since 
the normalized likelihood function is independent of density, 
the STY method provides two of the three LF parameters {a 
and M*) and is furthermore free of the problem of density in- 
homogeneities, provided that there is no correlation between 
the LF and density. However, for the same reason the STY 
method does not provide (jf, for which we need to resort to an 
independent method. 

4.2.3. Normalization for the STY-estimated LF 

To estimate the number density parameter 0*, we use the fol- 
lowing unbiased estimator for the mean number density « from 
Davis &Huchra( 1982) 



J s(z)w(z)dV ■ 



(15) 



Here, Wi(z) is a weighting function for galaxy / at redshift z, and 
s(z) is the selection function, which we define in redshift space 



rmm{L„iax{z)-:L2) 
Jmax(L\ ,Limn{z)) 



(j){L)dL 



rmdL 



(16) 



In this equation, Li and L2 are the minimum and maximum lu- 
minosities of the luminosity interval over which we would like 
to determine n, and L,„,„(z) and L„,ax{z) are the minimum and 
maximum luminosities observable at redshift z for given survey 
apparent-magnitude limits. 

The variance of this estimator is 



5n' = 



2 rP' j dV\dV2S\S2W\W2S,{r\Tr2) + n J dVsw 
(/ dVsw)^ ' 



(17) 



where the integral is done over the survey volume, and 
£,ifi,f2,z) is the two-point correlation function at redshift z- 

The optimal weighting function that minimizes the variance 
is roughly 

w = 1 / s{z) 



and we use this weight and equation ( |17[ ) to estimate the num- 
ber density and its error When the variance is minimized, the 
fractional error for the measurement of n is roughly (Davis & 
Huchra 1982) 



Sn 
n 



0.5 



(18) 



where 73 = 47r / r^£_{r)dr. For nearby galaxies, ^3 = lO'^/z'^Mpc^ 
(Lin et al. 1998). Note that this rough estimate based upon 
equation (jTsj) is accurate only when the depth of the survey 
volume in each dimension is much greater than the correlation 
scale (i.e., jc » rQ,y » ro, and z » ro). When the total volume 
is large but the depth of the volume in one or two dimensions 
is comparable to or less than the correlation scale (e.g., x < ro, 
and/or y < ro, as is our case here), equation (|l8|) overestimates 
the fractional error. 

To obtain a more accurate error estimate, we integrate equa- 
tion ( p^ ) numerically. This requires knowledge of the cluster- 
ing properties of E/SOs at the redshift of interest, which are 
not very well known. Nevertheless, as shown in Appendix B, 
we find it plausible to use an E/SO clustering evolution model 
with a spatial two-point correlation function that evolves as 
^(r) = {ro.E/so/ry'/(l+z)^-'^', with ro,£/so = 8 h-' Mpc (comov- 
ing coordinate), 7 = 1.8, and e = 0.8 and that cuts off at a scale 
r = 20/!-i Mpc. 

This clustering model gives 73 ~ 2 x 10^/i~^ Mpc^* from equa- 
tion (|l7|) at z = 0, while for the galaxy population as a whole at 
z = 0, the recent observed value is Jt, = 10'* /i"^ Mpc^* (Lin et 
al. 1996), with some previous estimates indicating a smaller 
/3 (e.g., 73=1700 from Davis & Peebles 1983). The fact that 
Ji ~ lO^/i""^ Mpc-^ from the analysis of the power spectrum and 
two-point correlation function of the Las Campanas redshift 
survey (LCRS) data indicates that the effect of the possible 130 
Mpc-scale structure (Landy et al. 1996) is negligible when 
estimating errors in the number density. Therefore the adopted 
cutoff at r = 20/1"' Mpc is well justified (see also Peebles 1994). 

With the assumed clustering model, we find that the frac- 
tional error of the number density of QS-E/SOs is about 20- 
25% at z ^ 0.8, and about 35% at z ^ 0.4, several times greater 
than the error estimates based on the Poisson statistics alone. 
These errors are comparable to the uncertainties and biases in 
the raw counts that were estimated in the previous section. This 
order of magnitude in the fractional error was also found by de 
Lapparent et al. (1989) for local galaxies. We have also tried 
increasing the cutoff value in the integral to r = 40 Mpc but 
find that the error estimate does not change significantly (at the 
10% level). 

Finally, the number density parameter, (/>*, is calculated as 



(19) 



Jl;'[(l)(M)/(i3*]dM 



where 



0(M)/0* 



is the partial LF derived from the STY method. 

Note that, since the value of (p* correlates with the estimate 
of M*, the uncertainty in M* introduces another source of error 
into 0*. We take this into account by adding the error con- 
tributed from M* to the error estimated above in quadrature. 
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4.3. Method! 

In this method, we parametrize the luminosity evolution as 



or equivalently. 



E{z) = Qz, 



M*(z)=M*(0)-E(z). 



(20) 



(21) 



The function E{z) is expressed in magnitude units, and thus 
galaxies get brighter as a function of redshift by E{z) mag. A 
similar formalism was used in the analysis of CNOC2 data by 
Lin et al. (1999). The linear form of equation ( pO| ) approx- 
imates the passively evolving stellar populations of the BC96 
models. 

For the number density evolution, we use the following pa- 
rameterization: 



\z)=4>*io)(i+zr 



(22) 



The value of m is claimed to lie between -1.5 and -1 for a 
CDM-dominated Einstein-de Sitter universe with hierarchical 
clustering (Baugh et al. 1996; Kauffmann et al. 1996). 

The parameters M* and 0* in equation (|l2|) are now replaced 
by M (z) (equation (|2l|)) and (f)*(z) (equation (p2|)), respec- 
tively, to provide the LF which incorporates the evolutionary 
change. 

In estimating the parameters of this LF (Q,m,M*{0),a and 
(f>*(Q)), we follow the procedure described in Lin et al. (1999). 
First, we use the STY method with defined in the same way 
as equation (|l4|) to estimate M*(0), Q, and a. The density evo- 
lution parameter (m) is then estimated by maximizing the like- 
lihood that each galaxy will be at its observed redshift, where 
the likelihood, L , is 



In L = E In (pi) + const ant, 



(23) 



and 



p, = p(Zi\MM,Q) = (1 +zd"'/ / (1 +2)"— ^2.4) 



dz 



Finally, (j>*iQ) is estimated analogously to the procedure in 
Section 4.2.3 except that we now take the density and the lumi- 
nosity evolution terms into consideration. To do this, we mod- 
ify equation (lla) as follows: 



/ s(z)w(z)dV 



(25) 



The absolute magnitudes in the selection function, s(z), are now 
calculated taking the luminosity evolution, E(z), into account. 
The error for (?!)*(0) is estimated using the method described in 
Section 4.2.3. 

4.4. The luminosity function ofE/SOs at z < 1 .2 

4.4. 1 . Results from Method 1 

Figure |T] shows the luminosity function of the Groth Strip 
QS-E/SOs in two redshift intervals, one at z = 0.05-0.6 (Zmed = 
0.45; figures on left with blue points and blue dashed lines) 
and the other at z = 0.6-1.2 (Zmed = 0.83; figures on right with 
red points and red sohd hues). Three different cosmologies are 



shown: top, Einstein-de Sitter universe; middle, an open uni- 
verse with rim = 0.2; and bottom, a flat universe with J7„, = 0.3 
and A = 0.7, which is currently favored (Im, Griffiths, & Rat- 
natunga 1997; Riess et al. 1998; Perlmutter et al. 1999; for 
a review, see Primack 2000). Also plotted is the local LF of 
E/SOs from Marzke et al. (1998) as a dotted line, and the local 
LF of E/SOs from Marinoni et al. (1999) as a dot-dashed line. 
Table 3 summarizes the LF parameters for the local LFs. Note 
that a ~ -1.0 in both cases, which justifies the use of a fixed a 
value for our fits. However, the M* values from the two works 
differ by about 0.5 magnitudes; the origin of this discrepancy 
may be due to differences in the magnitude systems, as dis- 
cussed in Marinoni et al. (1999). It is not clear which system 
is close to our magnitude system, thus we consider both values. 
We also notice that the value of (p* from Marinoni et al. (1999) 
is lower than the value from Markze et al. (1998) by nearly a 
factor of 2. This can be attributed to difference in classification 
schemes. Marinoni et al. (1999) select E/SOs as objects with 
-5 <T < -1.5, as opposed to Marzke et al. (1998), who take 
T <Q. 

The data points for the Groth Strip QS-E/SOs come from the 
1 /Vmax estimates, while the colored curves are the LF functions 
fit to the same points using the STY method. The bin sizes for 
the QS-E/SO LF points are determined so that each bin contains 
at least 3 galaxies, with the minimum bin size being 0.2 magni- 
tudes. Some LF points include as many as 15 E/SOs. 

To check for any strong biases in our sample (such as incom- 
pleteness), we calculate {V /Vmax) ■ The values of {V /Vmax) are 
found to be 0.56, 0.54, 0.53 (±0.03) for the Einstein-de Sitter, 
open, and A universes respectively. These values are close to 
those expected for passively evolving E/SOs 0.55), and thus 
imply both that the number of galaxies at high redshift is fairly 
complete (as found previously) and that the number density of 
E/SOs has not evolved dramatically. We return to the latter point 
after comparing 0* at z = with 0* at z ^ 1 quantitatively. 

Table 4 lists the best-fit LF parameters for the two distant red- 
shift intervals using Method 1 . Since the faint end of the LF is 
not well determined and since the number of E/SOs in our sam- 
ple is not large enough to provide a meaningful constraint on all 
three LF parameters simultaneously, we estimate M| by fixing 
a = -1 .0. This allows us to compare M| and (p* of Groth Strip 
QS-E/SOs with the local values in Table 3. Results are summa- 
rized in Table 5, which lists the shift in M| with respect to the 
local values, AB = M^iz = O)-M^(z). Here, the value of (p* 
from Marzke et al. (1998) in Table 3 has been multiplied by a 
rough correction factor of 0.7 in order to account for a possible 
difference in E/SO classification criteria (see below). 

A robust conclusion from Table 5 is that the bright end 
(L > L*) of the z ^ 0.8 LF is significantly brighter with re- 
spect to both local LFs, by about 1.1-1.6±0.2 mag relative to 
Marzke et al., and 0.7-1.2±0.2 mag, relative to Marinoni et al. 
This is most naturally explained by a luminosity brightening of 
all galaxies by these amounts from z = to 0.8. The precise 
amount depends slightly on the local samples and the assumed 
cosmology, but is largely independent of any biases or incom- 
pleteness in the actual counts. Moreover, the derived luminosity 
evolution is internally consistent, with AM| for the more dis- 
tant redshift bin being roughly twice as bright as in the nearer 
bin. Finally, the observed degree of brightening is much larger 
than uncertainties in our own magnitudes or any differences be- 
tween our own magnitude system and that of the local samples, 
which can amount to as much as 0.5 mag. 

We have calculated the expected magnitude brightening back 
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to z ~ 0.8 using passive-evolution models and varying both z/or 
and metallicity (CB96). The predicted amount of brightening 
at z ~ 0.8 from these models is comparable to the values listed 
in Table 5. In particular, for models involving short bursts of 
star formation followed by passive stellar aging, we find that 
Zfor ~ 1 -5 for AB ~ 1 .5 while z/«,- ^ 3 for AB ~ 1 .0. In the A 
universe, for which the brightening is largest, the bulk of stars 
in field E/SOs may have formed at redshifts as low as Zfor — 1 -5, 
especially if the local from Marzke et al. is used. For the 
open or Einstein-de Sitter universes, z/or > 2-3 is a better fit to 
the results. If alternative evolutionary models with extended 
star-formation times are used, the initial redshift of star for- 
mation can be pushed either earlier or later depending on the 
precise history of star formation. 

Related papers by our group discussing Groth Strip early- 
type galaxies estimate the amount of luminosity brightening by 
studying the fundamental plane of field E/SOs (Gebhardt et al. 
2000) and the size-luminosity relation of luminous bulges to 
z ~ 0.8 (Koo et al. 2000). The fundamental plane shows a 
brightening of 1 .4-1 .8 mag in the B band back to z ^ 0.8, while 
the size-luminosity relation of luminous bulges shows a similar 
result; both are consistent with the LF analysis here. Additional 
color information is used in these papers to further constrain 
spheroidal star-formation histories; the general conclusion is 
again that the bulk of star formation was quite early. 

The fits to (/)* in Table 5 can also be used to test for evolu- 
tion in the number density, assuming that the luminosity evolu- 
tion is well represented by the above change in M* and that the 
Groth Strip QS-E/SO sample counts identically the same kind 
of galaxies as the local samples. Our method tends to select 
E/SOs with T <-2 and may miss - 30% of E/SOs with T <Q 
according to the tests performed in section 3. Visual inspec- 
tion of red, non QS-E/SOs galaxies confirms that we may be 
missing roughly this many E/SOs. On the other hand, the local 
E/SO local sample of Marzke et al. (1998) is chosen to have 
r < 0, while Marinoni et al. is based on E/SOs with T < -1.5. 
This suggests that the Marzke et al. (1998) normalization may 
need to be multiplied by a factor of 0.7 to correct for classifi- 
cation mismatch. With this correction to Marzke et al.. Table 5 
shows that 0* at z ~ 0.8 is 1.6-2.0 times greater than the value 
at z = for the Einstein-de Sitter universe. This result is im- 
plausible as it indicates more E/SO galaxies at higher redshift 
than now, which is not predicted by any formation model. If it 
were well established, this trend might even suffice to rule out 
the Einstein-de Sitter cosmology, which in any case is highly 
disfavored by numerous other evidence (Bahcall et al. 1999; 
Primack 2000). However, given the ^ ±30% uncertainty in 
both the nearby and distant counts, the discrepancy must be 
considered marginal. The above increase in E/SOs at high red- 
shift disappears in the two low-il,„ cosmologies, for which the 
number density remains flat at all epochs within the errors. Al- 
though uncertainties are large, the data seem most compatible 
with a picture in which the great majority of E/SO galaxies were 
already in place by z ^ 1 and their numbers have been roughly 
constant since then. 

The above analysis using Method 1 relies on comparison to 
the local E/SO luminosity function, and thus admits numerous 
sources of systematic error, as we have noted. Another point 
not discussed so far is the highly correlated error between M* 
and (p* in the fitting procedure, i.e., when M* is estimated to be 
brighter, (/)* becomes smaller even for essentially similar data. 
Both of these problems can be minimized by using Method 2, 
which estimates the luminosity and number density evolution 



simultaneously from the Groth Strip QS-E/SOs alone. 

4.4.2. Results from Method 2 

In deriving the evolution parameters Q and m using Method 
2, we adopt two slightly different approaches. In the first ap- 
proach, we place no constraints on any of the LF parameters 
except for fixing a = -1.0. The derived M*(0) and 0*(O) are 
then compared to local values. In the second approach, we fix 
M*(0) to the value in Marzke et al. (1998) or in Marinoni et 
al. (1999). Fixing M*(0) helps reduce the error in Q and m, 
but the two local values of M*(0) differ by almost 0.5 mag. 
We consider both values and see what impact this has. Table 6 
summarizes the resultant LF parameters of E/SOs, with errors as 
indicated in parentheses. Figures ^^-^ show the significance 
of the measurements in contour plots of Q and m. As expected, 
larger values of Q are associated with more negative values of 
m. 

The first approach (solving for Mg(0) simultaneously) gives 
parameters for the local function that are consistent with the lo- 
cal data. The best-fit value forM^(O) is ~ -20.03 ± 0.6 mag for 
all three cosmologies, in almost perfect agreement with -20. 14 
mag as found by Marzke et al. (1998) but also consistent within 
the errors with the value -20.54 mag found by Marinoni et al. 
(1999). The estimates of also match well with the raw 

value of 0*(O) in Marinoni et al. (1999), or with the value from 
Marzke et al. (1998) corrected by the factor 0.70. 

The derived evolutionary parameters Q and m are also plau- 
sible and consistent with the results from Method 1 . The fitted 
Q values imply a luminosity brightening of 1.4-1.9 ±0.7 mag 
back to z = 1 in all models, showing definite evolution, as found 
also by Method 1 . The fitted values of m depend more strongly 
on cosmology, with m being ^ +0.5 for the Einstein-de Sitter 
universe (more objects at high redshift), versus -0.5 and -0.9 
for the open and A universes (fewer objects at high redshift). 
However, the error bars on m are ±0.7, so the data are in fact 
statistically consistent with no number density evolution in all 
three cosmologies; again, this agrees with Method 1. 

Tighter constraints on all parameters can be obtained by 
fixing Mg(0), as in the second approach. When we adopt 
MgiQ) = -20.14 from Marzke et al., we obtain almost the same 
values for the other parameters because the Marzke et al. value 
of Mg(0) is close to our own unconstrained value. However, the 
smaller error on Q now makes the detection of luminosity evo- 
lution much more secure (> 5-cr) under all three cosmologies. 
The error in the number density evolution also shrinks slightly, 
but the Einstein-de Sitter and open universe are still consistent 
with no evolution in number density, while the A universe fa- 
vors a slight decrease in past number densities at about the 2- 
(7 level (however, this result would no longer be statistically 
significant if our possible systematic undercounting of distant 
galaxies at the 30% level were included). 

When Mg(0) = -20.54 is adopted from Marinoni et al. 
( 1 999), we find a smaller amount of luminosity brightening due 
to the fact that this value of M|(0) for local E/SOs is brighter 
by 0.4 mag than that of Marzke et al. Nevertheless, the sig- 
nificance of the luminosity brightening at high z remains high. 
All values of m are increased by about 0.6 from the previous fit 
using Marzke et al., indicating more galaxies at high z than be- 
fore. Now, the Einstein-de Sitter universe shows an increase in 
number density at the lA-a level, while the other two universes 
are consistent with no number density evolution. 

To summarize. Methods 1 and 2 are basically consistent, in- 
dicating that the data are both internally consistent with them- 
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selves (low vs. high z) and externally consistent with the local 
samples. If the Einstein-de Sitter universe is excluded as un- 
likely, a brightening in Mg of 1.6-2.0 mag has been detected 
to z 1 . This result is robust and is largely independent of the 
count errors although it does depend on the relative definitions 
of the local and distant magnitude systems, which may differ 
by up to 0.5 mag. Evolution in the number density is much 
harder to determine given the errors; the data may favor a small 
decrease in number density at the level of a few tens of percent 
to z 1, but they are also consistent with no change. Sources 
of error in the number densities include raw and systematic er- 
rors in the distant counts, uncertainties in the local counts, and 
systematic classification differences between the local and dis- 
tant samples. Taken together, these inject errors of as much as 
a factor of two into the derived evolutionary trends in number 
density. 

4.5. Effect of changing redshift range 

As pointed out before, our results are potentially sensitive 
to a chosen redshift interval. We do not find bright E/SOs at 
z < 0.2 due to our bright magnitude cut, and this could poten- 
tially lead to an underestimate of the bright end of the nearby 
LP. There may be a genuine lack of bright E/SOs at z > 1 as 
well, in which case the bright end of LF (z ^ 0.8) may be 
also be affected. For these reasons, we have tried using dif- 
ferent redshift ranges to see if there is any strong change in 
the estimated evolutionary parameters. Since this is merely a 
sensitivity check, the cosmology is restricted to the open uni- 
verse for simplicity. Table 7 lists the results of the tests, for 
four different redshift Hmits: 0.2 < z < 1.2, 0.05 < z < 0.8, 
0.05 < z < 1.0, and 0.2 < z < 1.0, all using Method 2. Chang- 
ing Zmin changes the results by only a small amount {AQ = 0.1 
and Am ~ 0.2), in a direction that makes M* at lower redshifts 
brighter However, the effect of changing the upper redshift 
limit is more pronounced — reducing it makes the number den- 
sity parameter more positive by Am 1 or more with respect 
to the values listed in Table 6, while not changing Q very much. 
The sense of this confirms the earlier suspicion from Figure [I^ 
and[T^ that bright E/SOs may be missing beyond z > 1 .0 in our 
data. Possibilities for the cause of the apparent deficiency in- 
clude (i) the systematic underestimate of photometric redshifts 
at z > 0.8, and (ii) a genuine failure to detect bright E/SOs at 
z > 1 . Future analysis with near-lR data will be able to provide 
a firmer answer to hypothesis (i), as they enable better estimates 
of photometric redshifts at z > 0.8 by sampling the light above 
the 4000 break. Larger, deeper datasets would provide a check 
of hypothesis (ii). We also tried increasing the upper redshift 
limit, but find only negligible effects on the evolution parame- 
ter estimates. This suggests that observational selection effects 
act against detecting E/SOs at z > 1 .2. 

4.6. Monte Carlo simulation and Malmquist bias 

This section tests our results using a mock catalog created 
by Monte Carlo simulation for other possible biases due to ran- 
dom measurement errors. At z > 0.8, errors of E/SO photo- 
metric redshifts become greater, and the number of E/SOs in 
our sample declines as redshift increases. Hence, we expect 
to lose more galaxies below Zmax than adding more galaxies at 
above Zmax- This biases the number of high-z E/SOs to be un- 
derestimated. On the other hand, in a volume-limited sample, 
random, symmetric errors in redshift space make the measured 
redshifts to be always slighlty underestimated due to the fact 
that more galaxies scattered in from higher redshifts than from 



lower redshifts simply because of the volume effect. Again, 
this kind of Malmquist bias leads to an underestimate of the 
absolute fluxes. The random error in apparent magnitude cre- 
ates a bias in another direction. At around L* , as the num- 
ber of L < L* galaxies are greater than the number of L > L* 
galaxies, the random scattering in the flux would put more 
L> L* galaxies than there actually are, leading to overestimat- 
ing the bright end of the LF. Our LF analysis takes into account 
some of these biases by incorporating measurement errors in 
the LF estimates. However, our procedure is not perfect. In 
order to test any biases left in the LF measurements due to ran- 
dom measurement errors in apparent magnitudes and redshifts, 
we created mock catalogs of E/SO galaxies using the Monte- 
Carlo simulation procedure similar to that described in Im et 
al. (1995a). Into the simulated catalogs, we introduced ran- 
dom measurement errors. Then, we estimated the LF param- 
eters from mock catalogs, and compared the outputs against 
input LF parameters. By exploring a range of LF parameters 
and cosmologies, we find that AQ= Qmaput - Qinpui < 0. 1 , and 
< Am = moutput -I'Umput <; -0.35. This test shows that the 
random errors slightly overestimate the amount of luminosity 
and number density evolution. To be precise, the evolution pa- 
rameters listed in Table 5-7 would need to be corrected by this 
amount, but, since the amount of systematic bias is smaller than 
l-CT errors, we neglect this effect. 

5. DISCUSSION 

5.1. Implications for the merging rate since z = 1 

If E/SOs are formed via mergers of disk-dominated galaxies 
with similar masses (e.g., Kauffmann et al. 1993) and not de- 
stroyed or converted to spirals afterwards, the number density 
of E/SOs can only decrease as we look back in time. When the 
number density evolution is modeled as n(z) = n(0) x (1 +z)"' 
with n(0) being the number density of E/SOs at z = 0, semi- 
analytic models based upon an il„ = 1 CDM-dominated uni- 
verse predict m = -1.5 to -1.0, i.e, at least half of all mas- 
sive E/SOs today were formed via major mergers since z ~ 1 
(Baugh et al. 1996; Kauffmann et al. 1996). Our results 
for the Einstein-de Sitter universe are clearly inconsistent with 
such model predictions, excluding the value m = -1 (doubling 
in numbers since z = 1) at more than the 97-99.7% confidence 
level. 

A way to make the merging picture more consistent with an 
Einstein-de Sitter universe involves an alternative merging sce- 
nario wherein early-type galaxies simply increase their mass 
monotonically through minor mergers with other early-type 
galaxies or disk galaxies, leaving the overall number density of 
large E/SOs unchanged. This brightening would mimic AB evo- 
lution in the opposite direction to passive aging. Van Dokkum 
et al. (1999) find that close, bright pairs in one distant cluster 
mostly consist of red, early-type galaxies rather than late-type 
galaxies. If such a trend can be found to apply in low density 
environments, this would provide good observational support 
for the alternative merging scenario. 

When an open or a flat non-zero A universe are assumed, 
our count data admit the possibility of a moderate increase in 
number density from z = 1 to now. Taken literally, the best-fit 
values of m suggest an increase of a few tens of a percent in 
the numbers of E/SOs since z = 1 . Semi-analytical models tend 
to predict that number density evolution of early-type galaxies 
is weaker in open or A universes (e.g., Kauffmann & Chariot 
1998); our results are quite compatible with such predictions. 
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The little-or-no number density evolution seen for the mostly 
"red" QS-E/SOs further implies that the expected number of 
blue luminous E/SOs should be small. If brief, episodic bursts 
of star formation make E/SOs bluer (e.g., Chariot & Silk 1994; 
Trager et al. 2000 for evidence from local ellipticals), the in- 
terval of time for finding blue colors combined with a smooth, 
undisturbed morphology is only of order < 1 Gyr Assume that 
the duration of a burst is 5t, that the fraction of present-day 
E/SOs that underwent a starburst since z = 1 is //,„„;, that the 
average number of bursts per E/SO since z = 1 is iihurst, and that 
the number of bursts is roughly constant in time since z = 1 to 
now. The fraction of distant E/SOs that will appear blue (fbiue) 
is then 



f nnf.r ^^urst f burst . {5t/l Gyr) 
h,ue = 0.06 (^)(-5y)(,(^=i)/8Gy,y 



(26) 



where t{z= 1) is the look-back time from z = to z = 1. Our 
counts suggest that the fraction of blue E/SOs is not more than 
10% of all E/SOs, which is consistent with the small percentage 
predicted by equation (|6|). Note that this argument assumes 
that episodic star formation occurs randomly from z = 1 to z = 0; 
if it occurs predominantly at a particular redshift, we should see 
an increased number of blue E/SOs at around that redshift. Cur- 
rently, our data are not sufficient to test this. 

An alternative interpretation of the count data is that distant 
E/SOs consist at least partially of bulges of galaxies that will 
later accrete disks to become spirals. Comparison of the num- 
ber density of such objects at different redshifts would then be 
meaningless. However, this is quite unlikely, at least at the 
bright end of the LF, since the bulges of local spirals are not as 
luminous or massive as the brightest E/SOs. Recent works fur- 
thermore suggest that large disk galaxies were already largely 
in place at z ~ 1 (Vogt et al. 1997; Lilly et al. 1998; Simard et 
al. 1999) and that normal spiral galaxies out to z — 1 were as 
abundant as those at z = (Im et al. 1999; Brinchmann et al. 
1998; Driver et al. 1998). Wholesale transformation of galaxy 
types does not look likely at the present time but is a possibility 
that clearly must be studied further. 

Overall, the LF of the Groth Strip QS-E/SOs is quite con- 
sistent with the view that the majority of luminous E/SOs were 
already in place at z ~ 1 and that their luminosities have evolved 
smoothly and quiescently over time, with only a small number 
of significant star-formation bursts per galaxy since that epoch. 

5.2. Comparison with other studies 

The definition of our Q parameter for the luminosity evo- 
lution is identical to that in CNOC2 (Lin et al. 1999), which 
was defined for a sample at z < 0.55. Their sample was se- 
lected based on color rather than morphology, but our Q val- 
ues are in good agreement with theirs (2= 1 .58 - 2.00 ± 0.49), 
excepting the case where M^(0) in Marinoni et al. (1999) is 
assumed. For a morphologically selected distant field elliptical 
sample, Schade et al. (1999) measured a luminosity brightening 
of AMg = 0.97 ± 0.14 mag from z = to z — 1 in an Einstein- 
de Sitter universe. This again is in reasonably good agreement 
with our values. Using another morphologically selected field 
early-type galaxy sample, Im et al. (1996) reported a luminos- 
ity brightening of AM^ ~ (0.6 - 1.5) ± 0.5 mag back to z ~ 1 
for a flat universe with or without A, again consistent with the 
present results. As noted, the fundamental plane (FP) of high 
redshift early-type galaxies provides yet another independent 
estimate of luminosity evolution. For cluster early-types, the 



luminosity brightening is only about 0.75-1.0 magnitudes in 
rest-frame B back to z — 0.83 (van Dokkum et al. 1998). The 
brightening we find here for field early-types is larger than this 
but is consistent with the stronger luminosity evolution found 
for field E/SOs using the fundamental plane (Treu et al. 1999; 
Gebhardt et al. 2000). Thus all evidence seems to agree in im- 
plying significant luminosity brightening of field E/SOs back to 
z- 1. 

Results on number density evolution do not agree nearly as 
well. The CNOC2 group found rapid number density evo- 
lution proportional to lO" "*^-" with P = -1.07 ± 0.49 and P = 
-1.79 ±0.49 for the Einstein-de Sitter and open universes, re- 
spectively. If extrapolated, these values correspond to only 0.4 
and 0.2 times the present number density of galaxies at z ~ 1 . 
Kauffmann, Chariot, & White (KCW; 1996) likewise found that 
number density decreased strongly back in time, as (l+z)"'^ 
in an Einstein-de Sitter universe. These studies both disagree 
with our estimtes of low or non-existent evolution in Table 6. 
We have noted that both the CNOC2 and KCW samples were 
selected based purely on colors; we suspect this might be the 
cause of at least a part of the discrepancy, for the following 
reason. At low redshift (z < 0.4), the spread in colors is small 
for the different galaxy types, and it is easy for photometric 
errors to bump blue galaxies into a red-galaxy sample. The pre- 
cise choice of color boundary also matters greatly. Both effects 
could lead to an increase in in the apparent number density of 
red galaxies at low redshift. The CFRS sample used by KCW is 
furthermore known to be deficient in red objects beyond z > 0.8 
(Totani & Yoshii 1998; Im & Casertano 2000). 

Results from the morphologically selected samples of 
Schade et al. (1999) and Menanteau et al. (1999) appear su- 
perficially to agree better with the present work, but there is 
an important difference. Both works showed little evolution in 
claimed E/SO number density, but their samples include a sub- 
stantial number of blue E/SOs. Such blue objects are not present 
in our more tightly selected morphological sample, and our red 
objects by themselves are steady. 

The only previous data in genuine agreement with ours is 
the morphologically selected sample of Im et al. (1996). The 
work by Im et al. (1996) shows a luminosity brightening of or- 
der AB = 1 .73 mag out to z — 1 for a flat non-zero A universe. 
They also found little decline in number density to z = 1 , and 
concluded that more than 70% of z = E/SO galaxies seemed 
to be formed before z = 1 . These numbers match well with 
our results, although their sample may have been slightly more 
loosely selected than ours (they found 6.7 E/SOs per WFPC2 
field vs. 5.1 E/SOs per WFPC2 field in the present study). 

5.3. Is the density of E/SO galaxies in the Groth Strip typical? 

Since E/SO galaxies preferentially live in high-density en- 
vironments, our counts could be biased too high if the Groth 
Strip contains superclusters. We review here briefly some tests 
of this hypothesis in fields that are not adjacent to the Groth 
Strip. One such test looked at E/SOs in the Hubble Deep Field 
flanking fields (hereafter HDFF; Williams et al. 1996); no sig- 
nificant difference was found with regard to number density or 
other properties of E/SOs below z = 1. Another test checked 
the number of E/SOs in the GSS versus the number in the HST 
WFPC2 fields that were used for the studies of Im et al. (1996; 
1999). The mean surface density of E/SOs in the GSS is actually 
slightly smaller (a factor of 0.8) than the mean surface density 
of E/SOs in these 30 other HST WFPC2 fields, showing that 
the Groth Strip QS-E/SOs are not dominated by populations in 
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high-density regions. These two tests suggest that E/SOs in the 
Groth Strip can be considered as representative of field E/SOs. 

Also note that Fig. 8a shows several prominent peaks in the 
redshift distribution of GSS galaxies. A large fraction of E/SOs 
at z > 0.6 are associated with these peaks at z ~ 0.8 and z— 1.0. 
Thus, the derived amount of the evolution could be potentially 
biased due to the existence of these peaks. One test for this 
is to derive the evolution parameters using different redshift 
cuts. We estimate the evolution parameters adopting higher 
redshift cuts which excludes these peaks (e.g., 0.05 < z < 0.78 
and 0.05 < 0.95). We find that output m and Q values derived 
with these redshifts cuts are similar to the values listed in Ta- 
ble 7, suggesting that the peaky distribution itself does not af- 
fect the results significantly. We also point out that the rather 
wide adopted redshift interval 0.6 < z < 1 .2 for the LF (derived 
with method 1) averages out the peaky distribution. The amount 
of the evolution estimated with the "averaged-out" LF matches 
very well that derived with an alternative method (method 2), 
yet another indication that the peaky distribution is not affecting 
significantly our results. Certainly, a more decisive statement 
can be made for this issue by analyzing a larger set of data. 

5.4. Uncertainty in the local number density of E/SOs 

We have already mentioned various uncertainties regarding 
the normalization of the local LF, </>*. One issue is the pre- 
cise range of Hubble types in the various samples, which may 
account for the much of the difference of a factor of two be- 
tween Marinoni et al. (1999) and Marzke et al. (1998). An 
additional question is the normalization of the local LF for all 
types of galaxies together (see Marzke et al. 1998 and refer- 
ences therein), with the local universe appearing to be under- 
dense in some studies by as much as a factor of 1.8 (e.g., Elhs 
et al. 1996). Driver et al. (1996, 1998) fit the observed num- 
ber counts of E/SOs using the LF of Marzke et al. (1998), but 
with the local normalization boosted up by this factor. If such a 
renormalization is necessary, our conclusions regarding modest 
number density evolution would change drastically — the num- 
ber density of QS-E/SOs at 0.05 < z < 0.6 would be about two 
times smaller than the local value, meaning that about 50% of 
present-day E/SOs would had to have formed very recently! 

There is independent evidence that the high normalization 
adopted by Driver et al. (1996) is not necessary; they invoked 
it to fit the bright end of their E/SO number counts in Driver et 
al. (1995). However, Im et al. (2000) studied E/SO number 
counts using an expanded sample of 56 HST WFPC2 fields, 
which contained the 6 fields used by Driver et al. (1996, 1998) 
as a subset. Im et al. find that E/SOs in the 6 fields of Driver et 
al. are 1 .5 times more frequent than average, eliminating most 
of the discrepancy with the local normaUzation. We conclude 
that the normalization of the overall local LF by Marzke et al. 
is consistent and that the numbers of E/SOs in this function are 
reasonable after correction by the previously justified factor of 
0.7 to match our Hubble types. 

5.5. Colors 

Figure 8 indicates that the color-redshift relation of Groth 
Strip QS-E/SOs is in a reasonably good agreement with pre- 
dictions of passive luminosity evolution models. However, we 
find that the V -/ colors of Groth Strip QS-E/SOs at z > 0.8 are 
somewhat redder than the model predictions. Koo et al. (2000) 
and Gebhardt et al. (2000) analyze colors of early-type galaxies 
and luminous bulges and find that the rest-frame U —B colors 
(roughly, observed V-/) of these objects at z ~ 0.8 are 0.2-0.3 



mag redder than what passively evolving models predict for a 
luminosity evolution AZ? of ^ -1.8 mag or more. The problem 
is lessened by using a high-metallicity model, but it does not go 
away entirely. Other possible solutions include dust-extinction, 
more complicated star formation histories, and uncertainties in 
the stellar evolutionary models. For more discussion, see Geb- 
hardt et al. (2000) and Koo et al. (2000). 

6. CONCLUSIONS 

Using the residual parameter {R) and bulge-to-total light ra- 
tio (B/T), we successfully separate distant E/SOs from other 
types of galaxies. With this quantitative classification scheme, 
we identify 145 E/SOs at 16.5 < / < 22 in the Groth Strip that 
lie at z ;^ 1-2, termed "quantitatively selected" QS-ESOs. Spec- 
troscopic redshifts are available for 44 of these QS-E/SOs, and 
we find a very tight correlation between Zspec and V-I in the 
sense that their colors are the reddest among field galaxies at 
each redshift over the redshift range < z < 1.2. We use this 
tight correlation of color and Zspec to estimate redshifts for the 
remaining Groth Strip QS-E/SOs without spectroscopy and find 
that these photometric redshifts, Zphot, are accurate to 10% 
for Zphoi <, 1- 

Using the full sample of 145 Groth Strip QS-E/SOs, we con- 
struct their luminosity function at z — 0.8 (0.6 < z < 1.2) and 
at z 0.4 (0.05 < z < 0.6) for three cosmological models: an 
Einstein-de Sitter universe with = 1, an open universe with 
flm = 0.2, and a flat universe with il„, = 0.3 and A = 0.7. A ro- 
bust result is that rest-frame B magnitudes have brightened by 
AB ~ 1.1-1.9 mag since z = 1, with larger evolution taking 
place in open and A universes. This result is consistent with 
previous studies of the distant E/SO luminosity function and 
with our own studies of distant bulges and the field E/SO funda- 
mental plane (in preparation). Evolution in the number density 
of E/SOs is less well constrained, and pushing all errors to their 
maximum values yields uncertainties of a factor of two in num- 
ber density to z ^ 1 . For the open and A universes, the data are 
most consistent with roughly constant numbers of E/SOs back 
in time, perhaps favoring a moderate decline in numbers by a 
few tens of percent at z ^ 1 . For the Einstein-de Sitter universe, 
the data favor an increase, and any drop in galaxies by as much 
as a factor of two at z ^ 1 is strongly ruled out. 

The amount of luminosity evolution estimated from Mg at 
z ~ 0.8 implies that the major formation epoch of stars in E/SOs 
occurred at z/or ^ 2-3 in the Einstein-de Sitter and open uni- 
verse models, shifted down to as recently as z/or ^ 1 .5 in the A 
model. The large amount of evolution coupled with the mod- 
est change in number density is consistent with a picture in 
which the majority of luminous E/SOs galaxies in the field to- 
day (;> 70%) already existed at z ^ 1 and that they have not 
undergone dramatic evolution other than steady dimming of 
their stellar populations since then. If the major merging of 
disk galaxies is responsible for the formation of luminous field 
E/SOs, such a process must have happened predominantly be- 
fore z ^ 1. 
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Hubble Space Telescope, obtained at the Space Telescope Sci- 
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Fig. 1 . — Images of nearby E (top), Sbc (middle), and Im (bottom) galaxies and their residual images after subtracting the best-fit model image. The azimuthally 
averaged one-dimensional surface brightness profile along the major axis (points with en'ors) is also plotted and compared to the model-fit 1 -D profile for the bulge 
(dashed line), the disk (dotted line), and the sum of bulge plus disk (solid fine). Numbers below the residual images indicate the bulge-to-total light ratio (B/T) and 
the residual parameter (K). 
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Fig. 2. — The quantity R vs. B/T for 80 nearby galaxies from Frei et al. (1999). Red squares are for T-type, T < —3, green triangles are for T = —2, and blue 
stars are for —2<T< 0. Other types of galaxies are marked with crosses. The lower-right corner of the figure surrounded by a box represents the region where Es 
and SOs are found without much contamination from other types of galaxies. 
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Fig. 3. — The asymmetry parameter R of Frei et al. galaxies shrunk by different amounts, vs. input R values (r^ = 5 pix). Galaxies are binned by output r^. 
Output values of R are reasonably well correlated with input values down to = 2-3 px. 
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Fig. 4. — B/T of Frei et al. galaxies shrunk to different sizes, vs. input B/T values (r^; = 5 pix). Galaxies are biimed by output ru. Output values of B/T are 
reasonably well correlated with input values down to r^; = 2-3 px. 
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Fig. 5. — R vs. B/T for simulated Frei et al. galaxies shrunk to different sizes. Red squares are for galaxy types T < -3, green triangles are for — 3 < T < —2, 
and blue stars are for — 2 < T < 0. Other types of galaxies are marlced with blaclc crosses. Galaxies are binned by output r;,/. The lines at lower right comer of each 
plot represent a morphological cut at R = 0.08 and B/T = 0.4. Within this boundary, the sample of QS-E/SOs is somewhat contaminated by spiral interlopers for 
Oi/ ~ i-2 pix. However, this contamination is reduced in practice by reducing the R boundary for galaxies with very small apparent sizes, as described in the text. 
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Fig. 7. — Signal-to-noise {S/N) of objects in the Groth Strip. The S/N is calculated within a circle of radius = r«. The soUd line is a theoretical line for a galaxy 
with r« = 0.3 arcsec (3 pix). Objects that form the higher S/N sequence are from the deep GSS pointing with long expositfe time. 
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Fig. 8. — This figure shows the effect of each selection cut applied separately to isolate red-envelope E/SOs. V—I colors are plotted vs. spectroscopic redshifts for 
GSS galaxies with spectroscopic redshifts (the Zspec sample) and 16.5 < I < 22. Figure (a): All 262 objects with Zspec are plotted. The large, red circles represent 
objects with / < 20, the mid-sized, green circles show objects with 21 < / < 21.0, and the small, blue circles show objects with 21.0 < / < 22.0. Also plotted are 
the theoretical upper and lower color boundaries for passively evolving E-galaxy models (see text). Figure (b): Objects that survive the residual parameter (R) cut 
are plotted. Note that the R cut used is a function of the apparent half-light radius, as explained in the text. Figure (c): Objects that survive the bulge-to-total light 
ratio (B/T) cut are plotted. Figure (d): Final, "quantitatively selected (QS) E/SOs" are plotted, which consist of 44 objects that survive both the R and B/T cuts. 
The blue interlopers (filled circles) represent a failure of the selection technique. However, their number is small (15%), and their photometric redshifts based on 
V -I place them at low redshift, so they do not contaminate distant bins. Solid lines represent passively evolving solar-metallicity models with the formation times 
indicated. 
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Fig. 9. — This figure sliows how the QS-E/SO sample changes when the selection criteria are varied by small amounts. The total number of selected galaxies is 
shown in each panel; these are to be compared to the actual number of 44 galaxies selected in Figure 8. Figures (a) and (b) show how the sample changes when the 
B/T cut is changed, and Figures (c) and (d) show the effect of changing the R cut. Loosening the cuts mainly increases the number of interlopers, while tightening 
them mainly decreases the red target sample. Thus, the adopted cuts in Figure 8 are close to optimum. For details, see text. 
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Fig. 10. — Images of GSS QS-E/SOs in the Zspec sample, ordered by redshift. Numbers indicated in each panel are (from upper-left corner in clockwise direction): 
(1) GSS ID, (2) /, (3) V -I, (4) B/T, (5) residual parameter, R, (6) rest-frame B-band absolute magnitude (Mb), and (7) Zspcc- The 6 blue interlopers are plotted 
separately at the bottom of the panel. Displayed intensity values are roughly square-rooted ((Z)'"'^'^ to be precise) to show morphological features as clearly as 
possible (see text). The box size corresponds to 30 kpc in physical coordinates, assuming an open universe with f2 = 0.3 and h = 0.7. To compare with local samples, 
images of local E/SOs taken from Frei et al. (1996) are shown directly with the same pixel scaling. 
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Fig. 11. — Images of galaxies in the Zspec sample that are not classified as E/SOs (high R or low B/T), although their colors are as red as QS-E/SOs. Galaxy 
images are ordered by redshift. Numbers indicated in each panel are, from upper-left comer in clockwise direction: (1) GSS ID, (2) /, (3) V-I, (4) B/T, (5) R, (6) 
rest-frame B-band absolute magnitude (Ms), and (7) Zspec- 
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Spectroscopic z (z^p^^,) vs. Photometric z 




spec 



Fig. 12. — Spectroscopic redshifts (Zspec) vs. photometric redshifts (Zphot) for "red" QS-E/SOs. 



0.2 r 




Fig. 13. — The quantity Zspec vs. (Zphot-Zspec) tor the "red" QS-E/SOs of Figure 12. The solid hnes indicate the estimated rtns error, Szphot'> ^ function of Zspec- 
Note that dzphot increases when Zspec > 0.8 but stays within ^ 15% of Zspec out to Zspec = 1.2. 
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Fig. 14. — /814 vs. (V-/) of Groth Strip QS-E/SOs (squares) compared to all galaxies in GSS with Zspec- The thick squares are the "red" QS-E/SOs, thin squares 
are "blue interlopers," and crosses indicate remaining galaxies. Also plotted is a fiducial color cut that might be used for eliminating blue interlopers even without 
spectroscopic redshifts. We have not actually pruned the sample using this cut but have kept track of candidate blue interlopers using it. 
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Fig. 15. — Images of the 98 GSS E/SOs selected in the Zphot sample. Images are ordered by Zphot, excepting the last 11 objects, which are candidate blue 
interlopers. Numbers indicated in each panel are (from upper-left comer in clockwise direction): (1) GSS ID, (2) /, (3) (y-I), (4) B/T, (5) rest-frame B-band 
absolute magnitude (Ms), and (6) Zphot- 
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. Continued. 
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Fig. 16. — I,oi vs. z of GSS galaxies. QS-E/SOs ai'e mai'ked with squares, and blue interlopers (including candidates) are shown with triangles. Thick symbols 
represent QS-E/SOs in the spectroscopic sample, and thin symbols are QS-E/SOs in the photometiic redshift sample. Small crosses are non QS-E/SO galaxies with 
spectroscopic redshifts. Also plotted are lines of constant absolute magnitude (Mb =-21.14, -20.14, and -18.27). Solid lines include luminosity evolution of 1.7 X ^, 
while dashed lines assume no luminosity evolution. Note that Mb = —20. 14 corresponds to M* for the Marzke et al. local E/SO LF. 




Fig . 17. — Mb vs. z of QS-E/SOs. Meanings of symbols and lines are the same as in Fig.hd. Values of Mgs here are derived assuming no luminosity evolution. 
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DEEP E/SOs 
Es at z = 
SOs at z = 

Es + SOs at z = 0: 
Spirals at z = Oi 
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Fig. 18. — The ellipdcity distribution of Grotli Strip QS-E/SOs (heavy line). Also plotted (thin line) is a model distribution combining local cluster Es and SOs 
in which the fraction of Es is assumed to be 40%. The dashed line represents the contribution from Es. and the dotted line represents the contribution from SOs. 
The local sample of Es and SOs is taken from Dressier et al. (1980). The ellipticity distribution of spiral galaxies from Lambas et al. (1992) is also plotted as the 
dot-dashed line, for reference. 
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Fig. 19. — The ellipticity distribution of bright Groth Strip QS-E/SOs with Mb < -20 (heavy line). Also plotted is the ellipticity distribution combining local 
cluster Es and SOs in which the fraction of Es is assumed to be 80% (thin line). The dashed line represents the contribution from Es, and the dotted line represents 
the contribution from SOs. 
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Fig. 20. — The ellipticity distribution of faint Groth Strip QS-E/SOs with Mb > —20 (heavy line). Also plotted is the ellipticity distribution combining local 
cluster Es and SOs in which the fraction of Es is assumed to be 30% (thin line). The dashed Une represents the contribution from Es, and the dotted line represents 
the contribution from SOs. 



Table 1 

Adopted R cuts for selecting E/SOs from HST images 





m < 1 


l<rhi<2 


2 < rw < 3 


rm >3 


/<21 


<0.05 


<0.06 


<0.07 


<0.08 


21 </<21.5 


<0.05 


<0.06 


<0.06 


<0.06 


21.5 </<22 


<0.05 


<0.05 


<0.05 


<0.05 



Note. — r«s are in units of pixels. 
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Fig. 21. — The luminosity function of GSS E/SOs at two differenM^shift intervals for three different cosmologies. The LFs at 0.05 < z < 0.6 (zmei — 0.4) are 
indicated by blue points, dashed lines, and errors in the figures on the left. The LFs at 0.6 <z < 1.2 (zmea ~ 0-8) are shown by red points, dashed fines, and errors 
in the figures on the right. The red solid lines and the blue dashed lines are drawn using the LF parameters in Table 4, but we get virtually identical lines using the 
parameters in Table 5. To compare to z = 0, two local E/SO LFs are plotted — that of Marzke et al. (1998) with the dotted line, and that of Marinoni et al. (1999) with 
the dot-dashed line. 
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Mb,(0) varied during the fit 



Mb,(0) = -20.1, from Marzlte et al. 
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Mp,(0) = -20.5, from Marinoni ct al. 
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A=0.0 

h = 0.7 



12 3 4 

Q 

Fig. 22. — Error contours in m (number density evolution parameter) vs. Q (luminosity evolution parameter) for three different approaches for the Einstein-de 
Sitter universe. In panel A, Q and m are derived along with Mb» (0). In B and C, they are derived assuming fixed values of Mb» (0) from the local samples of Marzke 
et al. (1998) and Marinoni et al. (1999), respectively. The numbers indicated on each contour are Ax^ values. When projected onto the one-dimensional intervals 
either on m or Q, Ax^ = 1.00,2.71,4.00,6.63, and 9.00 correspond to confidence levels of 68.3%, 90%, 95.4%, 99%, and 99.73%. 
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M[3,(0) varied during the fit Mg,(0) = — 20. 1 , from Marzlce et al. 
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Fig. 23. — Same as Figure 22, for an open universe. 
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Table 2 

Quantitatively Selected E/SOs in the Groth Strip 



ID 



z code 



V 



(V-I) Mb B/T 



R 



r„(0."l) 



052_1555 0.9607 21.06 23.16 2.22 

052_6543 0.7555 1 20.31 22.15 1.89 

053_4446 0.5480 21.33 22.88 1.44 

053_6418 0.6750 1 20.98 22.66 1.81 

053_7537 0.6430 1 19.96 21.64 1.70 

054_2510 0.0810 1 16.59 17.53 0.93 

062_6859 0.9874 1 21.90 24.29 2.17 

064_3021 0.9965 21.13 23.09 2.06 

074_6044 0.9966 1 21.14 23.40 2.29 

082_6533 0.5029 20.57 22.04 1.46 

083_3771 1.2044 21.94 24.30 2.11 

083_6766 0.4322 1 21.00 22.34 1.32 

084_2525 0.8123 1 21.76 23.75 2.08 

092_1339 0.9031 1 21.40 22.70 1.23 

092_4957 0.2105 1 20.52 21.34 0.86 

093_2470 0.8110 1 19.69 21.53 2.09 

094_2660 0.9033 1 20.70 22.83 2.22 

094_2762 0.9326 1 21.31 23.19 2.00 

094_6234 0.8053 1 21.78 23.76 2.04 

102_2157 0.1445 21.77 22.72 0.99 

102_3649 0.5336 1 21.17 22.78 1.62 

102_5148 0.5665 20.79 22.36 1.47 

102_5177 0.6839 1 20.71 22.39 1.91 

102_5358 0.7190 19.87 21.69 1.92 

102_5558 0.7371 21.05 22.88 2.02 

102_6549 0.5752 19.10 20.72 1.69 

103_1115 0.4635 1 19.91 21.35 1.46 

103_4766 0.8118 1 21.10 23.14 2.09 

103_6061 0.3631 1 21.44 22.35 0.80 

103_7221 0.9013 1 20.79 22.83 2.31 

112_7644 0.5818 21.66 23.27 1.60 

113_3136 0.9447 21.97 24.06 1.91 

113_3311 0.8117 1 20.50 22.47 1.97 

114_1866 0.9328 21.50 23.57 2.05 

114_3760 0.7234 20.69 22.50 1.89 

122_7569 0.6896 21.18 22.94 1.76 

132_2723 0.6102 20.27 21.91 1.62 

133_1016 0.7572 20.88 22.74 2.00 

133_4463 0.6120 20.62 22.26 1.62 

133_6041 0.9674 20.82 22.93 2.13 

134_1633 0.3230 20.87 22.06 1.21 

134_4363 0.6287 21.97 23.66 1.60 

142_2752 0.3454 21.34 22.58 1.08 

142_3329 0.4891 20.09 21.54 1.55 

142_7077 0.0127 20.74 21.42 0.69 

142_7764 0.3188 21.79 23.01 1.04 

142_7871 0.2267 17.42 18.53 1.22 

143_2770 0.8045 21.39 23.32 2.06 

143_7957 0.4772 20.19 21.62 1.60 

144_3353 0.4492 20.94 22.33 1.36 

153_0711 0.3078 21.30 22.48 1.19 

153_2471 0.9873 21.68 23.81 2.05 

154 6532 0.7434 21.78 23.63 1.88 



-22.04 0.75 0.06 6.42 0.29 

-21.52 0.44 0.08 6.34 0.31 

-19.28 0.82 0.04 3.04 0.39 

-20.38 0.72 0.04 3.52 0.60 

-21.24 0.92 0.03 2.31 0.66 

-19.05 0.41 0.07 19.04 0.26 

-21.36 0.56 0.02 3.01 0.35 

-22.19 0.64 0.02 5.00 0.41 

-22.18 0.70 0.05 7.42 0.40 

-19.74 0.41 0.06 2.07 0.37 

-22.47 0.86 0.01 3.28 0.34 

-18.88 0.53 0.03 2.63 0.64 

-20.43 0.42 0.03 3.79 0.68 

-21.38 0.87 0.05 2.67 0.18 

-17.38 0.41 0.03 6.88 0.40 

-22.49 0.52 0.08 11.30 0.30 

-22.08 0.42 0.04 8.85 0.25 

-21.63 0.55 0.05 4.70 0.34 

-20.37 0.97 0.03 2.91 0.58 

-15.22 0.46 0.02 7.64 0.47 

-19.35 0.76 0.04 2.03 0.52 

-19.93 0.91 0.05 3.25 0.35 

-20.70 0.49 0.06 5.00 0.53 

-21.73 0.74 0.05 4.69 0.05 

-20.66 0.65 0.05 5.21 0.57 

-21.67 1.00 0.03 10.26 0.22 

-20.17 0.76 0.04 5.13 0.40 

-21.09 0.45 0.05 1.96 0.64 

-17.90 0.49 0.04 1.84 0.20 

-21.98 0.54 0.04 6.75 0.11 

-19.15 0.79 0.04 0.74 0.54 

-21.03 0.57 0.01 4.99 0.48 

-21.69 0.77 0.05 5.51 0.54 

-21.44 0.54 0.01 2.78 0.19 

-20.94 0.68 0.02 3.52 0.43 

-20.26 0.47 0.06 2.56 0.63 

-20.73 0.77 0.02 4.15 0.23 

-20.96 0.76 0.06 4.62 0.48 

-20.39 0.72 0.01 2.79 0.62 

-22.32 0.64 0.07 4.18 0.00 

-18.14 0.69 0.04 0.93 0.24 

-19.15 0.92 0.01 2.68 0.25 

-17.85 0.84 0.06 1.76 0.13 

-20.14 0.51 0.07 5.24 0.25 

-10.78 0.52 0.05 2.37 0.30 

-17.19 0.67 0.01 2.56 0.29 

-20.66 1.00 0.04 16.05 0.20 

-20.75 0.60 0.03 2.99 0.25 

-19.96 0.91 0.04 4.44 0.23 

-19.06 0.50 0.03 1.82 0.25 

-17.58 0.45 0.03 2.72 0.31 

-21.58 0.67 0.04 4.01 0.67 

-19.97 0.64 0.04 1.23 0.45 
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ID 



z code 



V 



(V-D 



Mb 



B/T R rw(0."l) 



162 


5349 


0.4804 


1 


20.11 


163_ 


5528 


1.0523 





21.79 


164 


1176 


0.0980 





20.70 


164_4063 


1.0060 





21.69 


164 


4638 


0.3198 





17.87 


164_ 


7956 


0.3868 





18.71 


172_ 


5049 


0.3564 


1 


21.18 


172 


7753 


0.6480 


1 


20.38 


173_ 


3911 


0.9798 





21.21 


173 


_4039 


0.5055 





18.60 


173 


4131 


0.2316 





21.29 


173 


.6766 


0.3954 





20.58 


174. 


,2027 


0.1279 





19.64 


174_ 


_7829 


0.8803 





21.89 


182 


4830 


0.6438 





21.92 


183 


1868 


0.8793 





20.41 


183 


2653 


0.7479 





19.74 


183 


_3056 


0.7748 





21.52 


183_ 


4478 


0.5865 





21.42 


192_ 


2330 


0.5727 


1 


20.65 


192. 


3330 


0.5759 


1 


20.61 


192 


_5343 


0.3532 





19.95 


193 


1227 


0.7992 


1 


21.18 


193. 


_1616 


0.3821 





18.99 


194 


6444 


0.2271 





19.72 


203 


.0833 


1.1566 





21.89 


203 


_1921 


0.5717 





21.50 


203 


_2022 


0.6708 





21.94 


203 


2622 


0.7505 





20.88 


203 


2634 


0.3660 


1 


17.83 


203 


3311 


0.7492 





21.84 


203 


3418 


0.6840 





20.96 


203 


5166 


0.3336 





20.51 


203 


5720 


0.6056 





20.73 


203 


7714 


0.6535 





21.42 


212_4836 


0.3522 


1 


19.35 


213. 


1764 


0.7904 





21.27 


213. 


_2362 


0.4439 





21.28 


213. 


6222 


0.1118 


1 


19.61 


213 


7564 


0.7105 





20.78 


214_ 


_2761 


0.7354 





19.51 


223_ 


_7276 


0.4338 





20.33 


224_ 


3977 


0.7824 





20.14 


224. 


4363 


0.4877 





19.50 


224_4413 


0.6695 





20.44 


232 


4429 


0.0463 





21.69 


233 


5568 


0.4242 





21.29 


233. 


.6431 


1.1180 


1 


21.08 


234_ 


4218 


0.5393 





20.98 


242_ 


5538 


0.4907 





21.55 


244 


4850 


0.3049 





19.98 


253 


1150 


0.5642 





20.10 


253 


2334 


0.5679 





21.33 


253 


_4345 


0.7043 





20.11 


253_ 


4452 


0.6119 





21.50 


254_ 


1671 


0.7492 





21.84 


262. 


2560 


0.6234 





21.13 


262 


7430 


0.6064 


1 


20.98 


263 


6262 


0.4098 





21.83 


263 


_6340 


0.3528 





21.30 


263_ 


6417 


0.4082 





21.67 


264. 


0931 


0.6780 





21.93 


264 


6053 


0.3955 





21.10 



21.61 
23.99 
21.52 
23.84 
19.12 
20.03 
22.25 
22.19 
23.33 
20.13 
22.35 
21.88 
20.49 
23.91 
23.63 
22.43 
21.61 
23.41 
23.03 
22.35 
22.25 
21.18 
23.19 
20.29 
20.74 
24.20 
23.09 
23.69 
22.73 
19.15 
23.70 
22.71 
21.71 
22.36 
23.13 
20.59 
23.18 
22.67 
20.46 
22.57 
21.37 
21.69 
22.05 
20.96 
22.17 
22.48 
22.65 
23.24 
22.51 
23.02 
21.13 
21.67 
22.91 
21.90 
23.15 
23.70 
22.79 
22.77 
23.19 
22.55 
23.02 
23.69 
22.41 



1.48 
1.97 
0.76 
2.12 
1.32 
1.39 
1.01 
1.75 
2.06 
1.63 
1.06 
1.30 
0.78 
2.04 
1.76 
1.96 
1.98 
1.92 
1.84 
1.73 
1.66 
1.21 
2.10 
1.28 
1.01 
2.25 
1.65 
1.80 
1.88 
1.37 
1.71 
1.85 
1.15 
1.51 
1.62 
1.25 
1.89 
1.39 
0.82 
1.79 
1.90 
1.51 
1.88 
1.50 
1.81 
0.77 
1.30 
2.07 
1.68 
1.46 
1.31 
1.51 
1.66 
1.83 
1.66 
1.99 
1.52 
1.76 
1.33 
1.20 
1.46 
1.90 
1.36 



-20.06 


0.73 


0.03 


3.65 


0.59 


-21.83 


0.65 


0.01 


2.52 


0.13 


-15.40 


0.46 


0.05 


4.10 


0.69 


-21.68 


0.82 


0.03 


3.88 


0.49 


-21.12 


0.59 


0.05 


12.23 


0.38 


-20.83 


0.63 


0.04 


7.74 


0.21 


-18.10 


0.73 


0.03 


1.87 


0.21 


-20.84 


0.86 


0.04 


2.55 


0.54 


-22.00 


0.59 


0.04 


4.78 


0.25 


-21.73 


0.58 


0.08 


7.01 


0.05 


-16.85 


0.69 


0.06 


7.88 


0.74 


-19.02 


0.46 


0.04 


2.38 


0.53 


-17.06 


0.55 


0.08 


3.02 


0.20 


-20.75 


0.71 


0.04 


2.88 


0.78 


-19.28 


0.84 


0.02 


1.77 


0.65 


-22.23 


0.54 


0.08 


6.17 


0.15 


-22.04 


0.53 


0.03 


8.98 


0.27 


-20.42 


0.65 


0.04 


1.96 


0.27 


-19.42 


0.68 


0.03 


3.87 


0.41 


-20.11 


0.63 


0.06 


6.82 


0.17 


-20.17 


0.41 


0.04 


2.29 


0.68 


-19.30 


0.90 


0.05 


5.10 


0.58 


-20.93 


0.73 


0.01 


3.64 


0.34 


-20.53 


0.82 


0.04 


5.77 


0.00 


-18.37 


0.88 


0.07 


2.15 


0.33 


-22.27 


0.70 


0.02 


4.95 


0.77 


-19.25 


0.86 


0.00 


2.23 


0.46 


-19.40 


0.62 


0.00 


1.29 


0.41 


-20.92 


0.93 


0.05 


3.02 


0.07 


-21.54 


0.83 


0.05 


11.64 


0.41 


-19.95 


0.51 


0.01 


3.14 


0.39 


-20.45 


0.87 


0.04 


3.10 


0.49 


-18.59 


0.94 


0.04 


1.70 


0.47 


-20.24 


0.70 


0.08 


4.73 


0.33 


-19.83 


0.50 


0.04 


6.20 


0.29 


-19.89 


0.51 


0.03 


4.31 


0.70 


-20.78 


0.44 


0.02 


2.09 


0.26 


-18.69 


0.76 


0.04 


1.26 


0.16 


-16.79 


0.41 


0.04 


6.38 


0.14 


-20.78 


0.83 


0.03 


2.71 


0.52 


-22.19 


0.56 


0.04 


7.06 


0.33 


-19.56 


0.73 


0.04 


5.98 


0.16 


-21.86 


0.71 


0.06 


2.36 


0.43 


-20.72 


0.56 


0.05 


5.04 


0.60 


-20.89 


0.64 


0.07 


4.58 


0.56 


-12.69 


0.71 


0.04 


6.54 


0.70 


-18.52 


0.44 


0.05 


2.52 


0.40 


-22.87 


0.79 


0.03 


5.60 


0.36 


-19.57 


0.48 


0.05 


4.42 


0.20 


-18.69 


0.61 


0.04 


1.91 


0.66 


-18.87 


0.53 


0.05 


3.65 


0.75 


-20.61 


0.60 


0.07 


4.25 


0.23 


-19.40 


0.77 


0.02 


1.39 


0.65 


-21.41 


0.53 


0.02 


5.64 


0.13 


-19.51 


0.96 


0.04 


4.00 


0.03 


-19.95 


0.61 


0.02 


3.67 


0.06 


-19.96 


0.98 


0.06 


3.22 


0.27 


-20.00 


0.83 


0.04 


2.12 


0.60 


-17.85 


0.82 


0.04 


2.03 


0.55 


-17.94 


0.46 


0.06 


3.77 


0.39 


-18.00 


0.57 


0.03 


1.27 


0.27 


-19.45 


0.57 


0.04 


1.41 


0.34 


-18.50 


0.66 


0.04 


4.82 


0.27 
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Table 2 — Continued 



ID 


z 


z code 


I 


V 


(V-I) 


Mb 


B/T 


R 


ruiO'.'l) 


e 


272_ 


_2255 


0.8842 





20.68 


22.70 


2.11 


-21.99 


0.87 


0.04 


4.37 


0.60 


272. 


.2871 


0.5063 





19.93 


21.41 


1.60 


-20.41 


0.51 


0.08 


7.44 


0.54 


272. 


^5241 


0.5053 





21.11 


22.59 


1.41 


-19.22 


0.52 


0.03 


1.55 


0.63 


273! 


!2671 


0.3945 





21.45 


22.77 


1.26 


-18.14 


0.59 


0.04 


2.77 


0.65 


273_ 


_5617 


0.4965 





20.59 


22.05 


1.48 


-19.68 


0.90 


0.03 


1.36 


0.70 


274_ 


_0837 


0.7059 





21.92 


23.72 


1.90 


-19.61 


0.87 


0.04 


1.85 


0.13 


274. 


.3875 


0.2826 


1 


19.49 


20.62 


1.12 


-19.15 


0.84 


0.04 


3.29 


0.67 


274_ 


_4341 


0.7137 





19.92 


21.73 


1.82 


-21.65 


0.50 


0.03 


4.56 


0.39 


274_ 


_5142 


0.6051 


1 


21.81 


23.39 


1.62 


-19.16 


0.49 


0.01 


2.62 


0.33 


274_ 


_5920 


0.8110 


1 


19.62 


21.66 


2.17 


-22.56 


0.48 


0.04 


10.74 


0.14 


282_ 


_5737 


0.7524 


1 


21.58 


23.45 


1.92 


-20.23 


0.75 


0.03 


3.17 


0.56 


283. 


.2254 


0.6504 


1 


20.14 


21.88 


1.77 


-21.10 


0.82 


0.06 


2.99 


0.72 


283_ 


_3250 


0.6509 


1 


19.87 


21.52 


1.87 


-21.37 


0.61 


0.05 


5.49 


0.11 


284_ 


_3854 


0.4495 


1 


19.56 


20.87 


1.38 


-20.44 


0.79 


0.04 


3.76 


0.53 


284_ 


5154 


0.2883 





18.93 


20.07 


1.18 


-19.77 


0.41 


0.06 


6.74 


0.32 


284_ 


6253 


0.2804 





19.75 


20.86 


1.14 


-18.87 


0.74 


0.05 


2.51 


0.10 


284. 


7275 


0.2661 





18.97 


20.07 


1.12 


-19.52 


0.53 


0.08 


4.68 


0.39 


292_ 


3076 


6540 


n 

yj 


?0 1 8 


21 89 


1 77 


—71 08 


88 

V7.00 


04 


7 86 

z.,o\j 


73 


292_ 


7235 


0.5034 





21.24 


22.72 


1.64 


-19.08 


0.80 


0.00 


5.67 


0.26 


294 


0718 


0.5156 


1 


19.37 


20.74 


1.52 


-21.03 


0.73 


0.03 


5.48 


0.56 


294. 


.2078 


0.9295 


1 


21.95 


23.39 


1.40 


-20.97 


0.59 


0.00 


5.38 


0.47 


294 


4544 


0.7484 





21.42 


23.27 


1.71 


-20.36 


0.46 


0.06 


4.26 


0.16 


302_ 


3631 


0.2569 





20.71 


21.79 


1.04 


-17.70 


0.92 


0.03 


2.21 


0.26 


312 


6405 


0.7965 





20.20 


22.13 


2.06 


-21.89 


0.66 


0.04 


5.31 


0.04 


313 


1250 


0.5148 


1 


20.14 


21.66 


1.58 


-20.25 


0.87 


0.05 


4.31 


0.37 


313 


3515 


0.7480 





21.46 


23.31 


1.78 


-20.32 


0.53 


0.03 


3.26 


0.43 


314 


2845 


0.1628 





21.94 


22.93 


1.05 


-15.35 


0.41 


0.04 


2.28 


0.25 



Note.— (1) ID: ID# of the GSS E/SOs; (2) z code:0 if Zphot, 1 if Zspec', (3) /: /-band total magnitude; 
(4) V: V-band total magnitude; (5) (V-I): V-I color; (6) Mb', rest frame absolute magnitude in B-band 
(see section 4.2. 1); (7) B/T: bulge-to-total light ratio measured in the /-band from our fits; (8) R: residual 
parameter measured in the /-band; (9) r^: half light radius measured in the /-band in units of pixels (one 
pixel=0."l); (10) e: eUipticity (see section 3.8) 



Table 3 

Parameters of Local E/SO LFs 



Sample 


a 


M*B 


(/)* (10-3 Mpc-3) 


Marzke et al. (1998) 


-1.00 ±0.09 


-20.14±0.10 


1.51 ±0.31 


Marinoni et al. (1999) 


-0.97 ±0.14 


-20.54 ±0.18 


0.84 ±0.20 



Note. — Parameters are adjusted assuming Hq = 70 km sec Mpc . 
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Table 4 

Parameters of LPs from Method 1 





= 


1,A = 




= 0.2, A = 0.0 


f2„ = 0.3, 


A = 0.7 


0.05 <z<0.6 


0.6<z< 1.2 


0.05 <z<0.6 


0.6 < I < 1.2 


0.05 < I < 0.6 


0.6 <! < 1.2 




-1.0 


-1.0 


-1.0 


-1.0 


-1.0 


-1.0 




-22.79 ±0.19 


-23.38 ±0.15 


-22,99 ±0.19 


-23 .72 ±0.1 5 


-23.17 ±0.19 


-23.92 ±0.15 




(1.35 ± 0.42) X 10"^ 


(1.70 ± 0.45) X 10"^ 


(0.99 ± 0.31) X 10" 


(1.01 ±0.29) X 10"^ 


(0.71 ±0.22) X 10"^ 


(0.68 ±0.19) X 10"' 


"S 


-20.62 ± 0.19 


-21.21 ±0.15 


-20.82 ±0.19 


-21.55 ±0.15 


-21.00 ±0.19 


-21.75 ±0.15 



Note. — Paramcicrs are calculated assuming Hq = 70 km sec Mpc . Errors include those caused by Poisson counting statistics and clustering statistics but not those caused by incompleteness or 
galaxy misclassifications. 



Table 5 

Luminosity and density evolution from Method 1 









1,A = 


fim ^ 0,2, 


A - 


iliii ^ 0,3, 


A ^ 0,7 




0.05 < z < 0,6 


0,6 < 1,2 


0,05 < z < 0,6 


0,6 < E < 1,2 


0,05 < E < 0,6 


0,6 < ; < 1,2 


AB(vs. Marzke) 
* (z) /(70% * (Marzke}} 
A5(vs. Marinoni) 
^ * (s}/ ^ * (Marinoni) 


0.48 ± 0,21 
1,28 ±0,48 
0,08 ±0,26 
1.61 ±0,63 


1,07 ± 0,18 
1,61 ± 0,54 
0,67 ± 0,23 
2,03 ± 0,72 


0,68 ±0,21 
0,94 ± 0,37 
0,28 ± 0,26 
1.18 ±0,46 


1,41 ±0,18 
0,96 ±0,34 
1,01 ±0,23 
1,21 ±0.45 


0,86 ± 0,21 
0,71 ± 0,26 
0,46 ± 0,26 
0.85 ± 0.33 


1,61 ±0,18 
0,64 ± 0,22 
1,21 ±0,23 
0,81 ±0,30 



Note. — The quantity '^*(0) from Marzke et al, (1998) has been multiplied by a correction factor of 0.70 to adjust for the fact that our galaxy selection criteria likely select 
fewer early-type galaxies than Marzke et al. See text. 



Table 6 

Parameters of LPs from Method 2 

nm = l,A = nm=0.2,A = f^m = 0.3, 



Solving for M*(0) 

Af*(0) -20,03 ±0,56 

Q 1.43 ±0,80 

III f) lf> ■- n !-N 

^* (1,10 = 0.44j X 10"- 



-20,05 ± 0,56 

1.82 ± 0.80 



-20,19 ±0,57 

1.89±0.81^ 



Fixed M* (0) = -20.14 from Marzke et al. (1998) 

Q 1,29 ± 0,23 1.70 ±0.23 



(1. 



0,64 ± 0,49 

II ±0. 



1-3 



-0,33 ± 0.49 
(1,02 ±0.35) X 10~- 



1,97 ± 0,23 
-0,95 ± 0.48 
(0.96 ±0,33) X 10~- 



Fixed M* (0) = -20.54 Iroiii Marinoni cl al, (1999) 

Q 0,76 ±0,23 1,17 ±0,23 

m 1,29 ± 0,52 0,33 ± 0,52 

^* (0.74 ± 0,27) X 10"^ (0.75 ± 0,27) X 10"- 



1,43 ± 0,23 
-0,30 ±0,51 
(0.71 ± 0,26) X 10"- 



Table 7 

Parameters of LPs for various z limits 



0.2 < z < 1.2 0.05 < z < 0.8 0.05 < z < 1.0 



Solving for M* (0) 

M*(0) -20,12 ±0,57 
Q 1,73 ±0,81 

m -0,68 ± 0,73 



1,99 ± 1,19 
0,96 ± 0,96 



-19,92 ± 0,57 -20,01 ± 0,72 -19,99 ± 0,58 
2,05 ±0,84 1,85 ±1,21 1,96 ±0,85 

0.16 ± 0,74 0,83 ± 1.07 0,00 ± 0,81 
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APPENDIX 

ERROR IN R FROM BACKGROUND NOISE 

Several effects contribute to the error estimates of the residual and asymmetry parameters (see Section 3, also Conselice et al. 
1999 and Wu 1999). Here, we quantify the error due to the background noise with the caveat that it is only a single component of 
the total error and therefore only a lower limit. 

We assume that the noise in all pixels is Gaussian, with true rms value a/,. This includes the object pixels, as background noise 
(from, for example, readout noise and sky photon statistics) is assumed to dominate everywhere. R is defined as 

R = Rt + Ra, (Al) 

where Rt and Ra are given in Eqs. (^ and (^. We further assume that both Ra and Rt are intrinsically zero, i.e., that the object is in- 
herently symmetric. It may be shown that only under this condition is the assumed background correction strictly valid. Furthermore, 
we are most interested in estimating errors for the case R is small, i.e., near zero. 

The (Ti, of the background is estimated locally because it may vary over the image. An error is introduced into Rj and Ra if (Jt is 
mis-estimated, as then the background correction, 

i:^-\B>j±Bf\, (A2) 

will be slightly wrong. However, with the noise being purely Gaussian locally, the only source of error in our estimate of ab comes 
from the fact that only a finite number of background pixels is available over which to estimate it. For this background estimate, we 
use the same number of pixels as the number of pixels under the object, namely, Npi_^. We calculate below the error in Eq. A2 due to 
the finite number of pixels in the background estimate. 

We assume that a mean background level has been subtracted, so that the distribution function of (B,y ± Bjj^^) is Gaussian with 

zero mean and a = Vlah- Given that Ra and Rj are both zero, it follows from the above assumptions that the distribution function of 
(Rij ± Rjj''''') is the same as that of (B,y ± B/j*"). Hence, there will be an identical error in the estimate of the quantity 

j:^-\Rij±Rf\. (A3) 

We can therefore add the two errors in quadrature to find the final total error in Ra and Rj due to Gaussian background noise. 

The error in the background correction, S(T,j\Bij zLBjj''^\), is simply Npi^ x 6b, where b is the estimated value of the background, 

{j\Bjj±Bjj'^\), and Sb is its error. The distribution function of j\Bij±B]j^^\ is not Gaussian, but the error of its estimated mean can 
nevertheless be calculated using the central limit theorem, provided that the number of samples (pixels) is large enough. The central 
limit theorem says that the computed mean using samples from a distribution F(x) approaches < jc > with an error ax/VN, where 
ax is the variance of x. The minimum for this to hold with some accuracy is typically = 100, which is generally true in our case 
since Npix > 100. The error in the background correction can then be shown to be 



S(j:-\Bij±Bjf'\)=Jj:iS\Bij\f, (A4) 



2 2 2 

a:=<x >-<x> . 



where (S\Bij\)^ is the variance of In deriving Eq. (A4) we have taken into account that the left side sums over all pairs of pixels 
twice. 

The formula for the variance is 
Applying this to the variance (S\Bij\)^ we have 
The first term in Eq. (A5) is 



(S\B,j\y = {\B,j\')-{\B,\y. (A5) 




(A6) 

Jo 

and the second term in Eq. (A6) is 

= / xexp dx/ \/ jl ab = \\ — ai,. (A7) 

Jo V TT 

Thus, the rms value of IB, , I is 



S\Bij\ = J^c7h. (A8) 
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Using Eq. (A4) and Eq. (A7), the uncertainty in the background correction is then, 



Si^i:\Bij±Bf\) = ./N~5\Bij\ = \j^Npi, X a^. (A9) 



Adding the identical term in quadrature for the residual sum, we obtain the total error in Rj and 



A, 



5RTorA=^ . (AlO) 

We ignore the error in the denominator — the sum in the denominator is a sum of total intensities whereas both sums in the numerator 
are sums of residuals, and they are furthermore subtracted, so it is fair to assume that the fractional error in the denominator is small 
and can be ignored. 

Since the signal-to-noise of an object image is defined as 

we have 

5RT„rA^Q-i5{SlN)-\ (All) 

However, R is the sum of Ra+Rt, and the pixels used are common to both Rt and Ra- Hence it is reasonable to assume that the 
errors 5Ra and SRj add arithmetically rather than in quadrature, and we finally get 

SRc^l.liS/Nyy (A12) 
One can also derive a similar relation for the asymmetry parameter. A, defined as 

mj-^}r\-nB,j-Bf\ 



where X/j is the flux or a flux related quantity at (/, j). In that case, we get an almost identical formula to Eq. A12, namely, 

SA = lJ{S/NrK (A13) 

SIMULATED IMAGES OF FREI ET AL. GALAXIES 

As we explained in Section 3, we rescaled and added noise to the images of local galaxies from Frei et al.(1996), in order to test 
our quantitative scheme of morphological classification when the image quality is close to the GSS data. The size-magnitude relation 



of Groth Strip QS-E/SOs is plotted in Figure B25 together with the size-magnitude relation of other types of galaxies. Like nearby 
galaxies, E/SOs tend to be more compact than other types of galaxies at a given magnitude since E/SOs have centrally concentrated 
light profiles. The r/,; ~ 5 pixels corresponds to rhi,med of / < 21 E/SOs, and r/,/ ~ 3 pixels is roughly equal to rhij„eci of 21 < / < 22 
E/SOs. Hence, images are rebinned and box-averaged so that the simulated galaxies have r/,/ ~ 5,3,2, 1.5, and 1 pixels. Noise is 
added so that the simulated galaxy images have S/N ~ 30-70, comparable to GSS galaxies at 20 < / < 22 (see Figure We 
also add an additional background area to the simulated image to help determining the background correction for R more accurately, 
since some of the Frei et al. galaxies almost fill up the original image, making it difficult to estimate the background noise. We 
have not convolved the images with the WFPC2 PSF; however, this does not affect significantly our results since the Wide Field 



Camera portion of the WFPC2 undersamples the PSF. Simulated images of representative galaxies are shown in Figure B26 When 
fhi ~ 3 pixels, the morphological details of late-type spiral galaxies are still somewhat visible, and R gives us a good measure of 
the complexity in their morphology. When r/,; ~ 2 - 3 pixels, it becomes challenging to classify galaxies by eye: Non-smooth, and 
asymmetric features of galaxy morphology are visible for late-type galaxies, but more detailed morphological features such as spiral 
arms are almost completely wiped out. For those objects, R still provides a reasonable measure of morphological detail from the 
non-smoothness and asymmetry in galaxy SB profile. When r/,/ ~ 1 pixels, it is almost impossible to classify galaxies by eye, but/? 
can still be used to identify the late-type galaxies with large R. 



CLUSTERING PROPERTIES OF E/SOS OUT TO Z ~ 1 

Studies of local galaxies show that bright E/SOs are more strongly clustered than other types of galaxies, and that E/SOs Uve 
preferentially in high density environments (Willmer, da Costa, & Pellegrini 1998; Loveday et al. 1995). However, there is no study 
which estimates the clustering of morphologically selected E/SOs at high redshift. In order to estimate the clustering of E/SOs, we 
can take three approaches: (1) use the clustering properties of bulge-dominated galaxies (e.g., Neuschaefer et al. 1997), (2) use the 
clustering properties of red galaxies as estimated from other redshift surveys (e.g., LeFevre et al. 1996), or (3) estimate the clustering 
properties of the Groth Strip QS-E/SOs directly from our GSS data. The first two approaches may suffer from contamination: in the 
first from late-type galaxies with high B/T (> 0.3), and in the second from early-type spirals as well as dusty late-type galaxies. 
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According to Neuschaefer et al. (1997), the amplitude of the two-point angular correlation function of bulge dominated galaxies 
(A„ in uj :^ A„ 9~^ '^) normalized at 1 ' is about 0. 1 at 1 8 < / < 22. Note that this value is several times greater than the value of A„ 
for the rest of the population in Neuschaefer et al. (1997), consistent with the clustering of E/SOs at z = relative to the other types 
of galaxies. Neuschaefer et al. (1997) also find that the A„ of MDS bulge dominated galaxies is well fit by a clustering model with 
a clustering scale of ro = 5.5 Mpc fixed in physical (or proper) coordinates. If we model the evolution of the spatial correlation 
function as £_(r) = {ro/ry /( \ +zy'~''^'^ with respect to the comoving coordinate system, the above case corresponds to e = 0. With 
this model, ro(z = 0.75) = 4 Ir^ Mpc in comoving coordinates; however, since studies of local E/SOs suggest that ro.£ /so ~ 6 - 8 
Mpc (Loveday et al. 1995; Guzzo et al. 1997; Willmer et al. 1998), it seems more reasonable to model the clustering evolution as 
^oiz) = ''o,£/so/(l +z)*^~'>'"^^^/t with e = 0.8. With this model, we also get the observed ampUtude of the two-point angular correlation 
function, A„ ~ 0. 1 . Furthermore, this model is consistent with the expectation of the linear growth of the density perturbation. 

As for the second approach, LeFevre et al. (1996) estimate that ro = 2.5 Mpc for red galaxies in the CFRS. This value appears 
to be much lower than the value quoted for bulge-dominated galaxies. In LeFevre et al. (1996), red galaxies are selected as galaxies 
with colors redder than those of Sbc type SED; it might be that more weakly clustered galaxies are introduced into the red galaxy 
sample by this somewhat less conservative color selection. 

The third approach seems quite attractive, since the clustering property is internally determined from the data we are studying. 
The variance of the number density of galaxies for a given survey geometry for which V '/^ < ro is 

Sn^ = n^i{0 + l/n) (CI) 
Here < ^ > is the average of the two-point correlation function over the volume, and can be approximated for a slab geometry as 

(0^3.9:^. (C2) 

To obtain the mean and the variance of the number density of GSS spheroids for a small volume element V, we divide our GSS 
volume into 10 equal volumes from z = 0.5 to z = 1 in the radial direction. From these volume elements, we find the Sn^ ~ 340, and 
n = 8.2. The clustering evolution as a function of redshift is ignored here for this order of magnitude calculation. With these values, 
we get ro 4h^^ Mpc (in comoving scale) at z 0.75, consistent with the value preferred by the study of the angular correlation 
function of bulges (Neuschaefer et al. 1997). Note that these values are estimated from the Groth Strip QS-E/SOs with spectroscopic 
redshifts by applying a selection function similar to that in Koo et al.(2000). The photometric redshift sample can not be used for this 
particular estimate of (^) since the uncertainty in the photometric redshift {6z ^ 0.1) is big enough to smooth out the true clustering 
property along the radial direction. 

The excellent agreement between the first and the third methods demonstrates that it is reasonable to assume the clustering of 
E/SOs as ^(r,z) = (ro/r)T/(l +z)^~'''^'^ where we adopt ro = 8/;"' Mpc, 7 = 1.8 and e = 0.8. Therefore, we adopt this clustering model 
for the Groth Strip QS-E/SOs to estimate the error in their number density. 
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Fig. B26. — Images of 1 6 representative Frei et al. galaxies, simulated to have r;,; ~ 5,3.2, and 1 pixels and S/N ~30-70. Images are ordered by morphological 
type (T-type). Numbers indicated in each panel are (from upper left corner in clockwise direction): (1) galaxy name (NGC number), (2) T-type, (3) output rf,i 
(pixels), (4) R, and (5) B/T. 
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